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PREFACE 


The  values  of  the  field  strength  given  in  the  section  on 
numerical  comparison  (5*1  and  5*2)  of  the  original  publication 
of  this  report  were  found  to  be  in  error.  Dr.  Paul  S.  Swan* 
first  suggested  the  possibility  of  an  error  and  Dr.  James  R. 

Wait**  showed  that  such  an  error  existed.  Subsequently  two  in¬ 
dependent  evaluations  of  these  numerical  values  were  made.  These 
were  in  close  agreement.  The  results  of  these  evaluations  are 
published  in  this  "corrected  edition"  of  Technical  Report  EE-44. 
Using  these  revised  field  strengths,  comparison  of  the  relative 
efficiencies  for  the  submerged  horizontal  and  vertical  dipoles 
shows  that  the  horizontal  submerged  dipole  is  more  effective 
than  the  vertical  submerged  dipole  in  exciting  a  vertically- 
polarized  field  in  the  air. 

For  a  receiving  site  located  1000  km  from  a  transmitting 
antenna  with  100  kw  of  input  power  at  1  cps,  a  signal-to-noise 
ratio  of  approximately  -15  db  is  predicted  by  the  mode  solution. 

At  closer  distances  the  branch-line  solution  indicates  that  tbs 
signal  is  considerably  stronger,  thus  suggesting  the  possibility 
of  practical  communications  at  distances  less  than  1000  km. 

This  "corrected  edition"  of  the  report  supersedes  the 
original  report  (Technical  Report  EE -44),  dated  February  1961, 
which  should  be  destroyed  and  replaced  by  this  "corrected  edition. 

3f 

Operations  Evaluation  Group,  MIT. 

National  Bureau  of  Standards,  Boulder,  Colorado. 


1.0  INTRODUCTION 


The  radio  wave  propagation  problem  as  it  exists  in  con¬ 
nection  with  communications  involving  submerged  submarines  is 
necessarily  one  which  concerns  low  frequencies.  This  conclusion 
is  the  immediate  result  of  considering  the  attenuation  of  an 
electromagnetic  wave  as  it  propagates  through  a  conducting  medium, 
and  noting  that  such  attenuation  is  proportional  to  a  factor 
exp  ( -cJ~E  x)  where  c  is  a  constant  for  the  medium,  f  is  fre¬ 
quency  and  x  is  distance.  When  x  is  as  much  as  a  few  meters 
and  the  medium  concerned  is  sea  water,  this  factor  automatically 
excludes,  for  practical  purposes,  any  frequencies  above  the  VLF 
and  ELF  categories.  It  is  customary  to  define  the  ELF  range  of 
frequencies  as  extending  roughly  from  1  to  3000  cps"*"  and  VLF 
from  3,000  to  30,000  cps . ^ 

In  the  initial  part  of  the  following  investigation,  ex¬ 
pressions  will  be  obtained  for  the  Hertz  potential  in  three  re¬ 
gions  of  arbitrary  conductive  and  dielectric  properties.  The 
regions  are  separated  from  each  other  by  two  infinite  parallel 
planes.  In  directions  normal  to  these  planes,  the  regions 
above  and  below  the  intermediate  layer  are  considered  to  be  of 

1J.  R.  Wait  and  Nancy  F.  Carter,  Field  Strength  Calcula- 
tions  for  ELF  Radio  Waves,  NBS  Tech.  Note  No.  Boulder 

Laboratories,  March  I960. 

2 

J.  R.  Wait,  A  Survey  and  Bibliography  of  Recent  Research 
in  the  Propagation  of  VLF  Radio  Waves,  NBS  Tech.  Note  No. 

Boulder  Laboratories,  May  1950. 
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semi-inf inite  extent.  The  regions  are  individually  homogeneous . 

An  electric  dipole  is  located  in  one  of  the  semi-infinite  media, 
which  will  be  called  medium  1.  The  intermediate  layer  is  called 
medium  2,  and  the  other  medium  J> . 

Following  this  initial  formulation,  specialization  will  be 
made  to  the  case  where  medium  1  is  sea  water,  medium  2  is  air, 
and  medium  3  is  ionosphere.  Under  these  conditions  the  iono¬ 
sphere  will  form  the  upper  boundary  of  a  parallel  plane  wave¬ 
guide  having  for  its  lower  boundary  the  surface  of  the  sea. 

The  theory  that  is  developed  is  valid,  up  to  a  certain  point, 
for  antennas  buried  in  earth.  The  point  at  which  a  departure  is 
made  to  the  specific  case  of  sea  water  is  in  the  evaluation  of 
the  branch  line  integrals,  where  certain  approximations  are  made 
based  upon  the  magnitudes  of  the  parameters ,  It  would  not  be 
difficult  to  perform  the  different  evaluations  necessary  for 
antennas  in  earth,  but  this  has  not  been  undertaken  in  this  paper. 

The  choice  of  a  parallel  plane  geometry  is  dictated  by 
reasons  of  simplicity.  There  is  no  intrinsic  reason  why  equations, 
however  involved,  cannot  be  worked  out  for  a  spherical  model. 

There  are,  however,  two  points  which  argue  against  developing 
the  more  formidable  case  first.  One  is  the  practical  difficulty 
of  extracting  numerical  information  from  the  results,  which 
would  involve  the  poorly  tabulated  spherical  Hankel  functions 
of  complex  order  and  argument .  The  other  is  the  fact  that  the 
parallel  plane  development  will  give  a  reasonably  good  estimate 
out  to  distances  from  the  antenna  on  the  order  of  the  earth's 
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radius,  and  this  is  an  ample  range  from  which  to  decide  whether 
a  spherical  model  approach,  required  for  the  greater  ranges,  is 
practically  justified.^ 

In  the  following,  both  vertical  and  horizontal  submerged 
electric  dipoles  are  considered.  In  each  case  a  "mode"  solution 
is  obtained,  corresponding  to  TEM-type  propagation  in  the  parallel 
plane  waveguide  for  which  the  boundaries  are  water  and  ionosphere . 
Restriction  of  the  frequency  range  to  1-1000  cps  avoids  the 
necessity  of  discussing  the  higher  order  modes,  which  are  then 
so  heavily  damped  as  to  be  of  no  practical  significance. 

Extension  of  the  results  to  higher  frequencies  is  nonethe¬ 
less  readily  carried  out,  provided  only  that  the  complex  charact¬ 
eristic  numbers  (occurring,  in  our  development,  as  poles  in  the 
complex  plane)  can  be  found  for  the  modes  in  which  appreciable 
transmission  can  occur. 

Investigators  at  the  Boulder  Laboratories  of  the  National 
Bureau  of  Standards  have,  apparently  been  successful  in  finding 
the  important  characteristic  numbers  in  connection  with  standard 
VLF  and  ELF  mode  propagation  theory.^  Since  the  modes  are  the 
same  no  matter  how  they  are  excited,  the  pole  locations  necessary 
for  consideration  of  higher  frequencies  are  therefore  available 
and  immediately  applicable  to  the  results  of  this  paper, 

^Howe,  H.  H.  and  Wait,  J.  R.,  "Mode  Calculations  for  VLF 
Ionospheric  Propagation",  Proc,  Symposium  on  Propagation  of 
VLF  Radio  Waves,  V.  5,  paper  no.  30,  Boulder,  Colorado,  1957 . 
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In  addition  to  fields  from  the  mode  solutions,  there  are 
field  contributions  due  to  the  integration  around  two  branch  cuts 
in  the  complex  plane.  The  relative  contributions  from  these 
branch  cuts  can  be  shown  to  diminish  rapidly  with  increasing 
distance  from  the  real  axis  of  this  complex  plane,  this  distance 
being  proportional  to  the  square  root  of  conductivity  times  fre¬ 
quency.  In  the  pages  that  follow,  we  find  that  the  integral  for 
one  of  these  cuts--that  connected  with  the  propagation  constant 
for  a  plane  wave  in  sea  water--is  practically  negligible  at 
large  distances  from  the  antenna.  Physically  this  means  that 
propagation  directly  through  the  water  from  transmitting  to  re¬ 
ceiving  antenna  can  be  disregarded. 

The  contribution  from  the  branch  cut  associated  with  the 
propagation  constant  for  a  plane  wave  in  the  ionosphere  is  not 
important  at  distances  out  to  around  1000  kilometers  for  the 
lower  part  of  the  ELF  range.  To  a  first  approximation  this  con¬ 
tribution  behaves  asymptotically  as  e~a^/p^,  where  p  is  dis¬ 
tance  in  meters.  The  factor  is  contained  in  the  quantity  a, 

which  is  the  imaginary  part  of  the  propagation  constant  k^,  of 
the  ionosphere.  Physically,  this  solution  corresponds  to  pro¬ 
pagation  along  the  air-ionosphere  boundary,  but  in  the  ionosphere. 

It  is  possible  to  draw  some  comparisons  between  the  results 
contained  herein,  and  those  which  have  been  obtained  by  other 
investigators  for  the  situation  when  the  antennas  are  located 
in  the  atmosphere.  In  the  latter  case  the  chief  interest  has 
been  in  VLF  propagation,  for  which  the  branch  line  integrals 
are  quite  rapidly  attenuated  with  distance  from  the  source, 
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and  therefore  can  be  neglected  in  comparison  with  the  mode 
solutions.  In  some  cases,  however,  the  VLF  mode  solutions  have 
been  extended  to  ELF  without,  apparently,  a  re-evaluation  of 
the  branch  line  contributions.  The  writer  has  given  considerable 
attention  to  this  matter,  but  the  results  are  necessarily  intended 
to  serve  only  as  estimates  when  the  product  of  distance  times 
square  root  of  frequency  is  sufficiently  great.  It  is  believed 
that  these  estimates-- in  spite  of  the  uncertainty  associated 
with  them- -have  served  the  very  useful  purpose  of  showing  that 
the  branch  line  solution--at  least  in  the  case  of  the  horizontal 
submerged  dipole- -is  behaving  in  such  a  way  as  to  merit  close 
attention  in  the  low  frequency  part  of  the  ELF  range,  at  dis¬ 
tances  less  than  1000  km. 

This  possible  importance  of  the  branch  line  solutions 
justifies  a  more  thorough  investigation  of  their  magnitudes. 

A  different  approach  is  then  required,  and  it  appears  that  a 
possible  line  of  attack  is  to  place  the  branch  cut  from  k^  at 
an  angle  of  (in  the  direction  of  k^  —  see  Figure  2.2). 

Then  tabulated  values  of  the  Bessel  functions  of  argument 
(e  iir/4x)  may  be  used  to  effect  a  numerical  integration.  This 
work  would  perhaps  be  best  performed  by  means  of  an  automatic 
computer.  It  will  be  found  later  (Section  6.0)  that  there  may 
be  some  hope  of  useful  communications  between  submerged  antennas 
at  1000  km.  This  preliminary  conclusion  is  based,  however,  on 
some  rather  optimistic  data  which  has  been  published  for  natural 
noise  in  this  frequency  range.  The  mora  exact  evaluation  of  the 
branch  line  integrals  should  therefore  perhaps  be  confined,  at 
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least  initially,  to  distances  of  1000  km  or  less,  which  might 
constitute  a  less  borderline  situation. 

It  has  been  found  that  for  a  submerged  short-circuited 
horizontal  coaxial  antenna  540  meters  long  and  a  transmitter 
power  of  100  kilowatts,  the  signal  power  at  1000  km  distance 
can  be  expected  to  be  about  13  db  below  noise  power  for  1  cps 
bandwidth  and  a  center  frequency  of  1  cps  excluding  the  exponen¬ 
tial  depth  factor.  This  bandwidth  is  quite  large  at  this  fre¬ 
quency;  hence  improvement  is  possible,  depending  on  how  much 
bandwidth  is  required  for  the  information  being  communicated . 

At  a  distance  of  100  km,  the  signal  power  would  be  roughly  40  db 
greater  due  to  the  mode  solution;  an  additional  contribution 
might  come  from  the  branch  line  integral  as  well. 

A  few  words  should  be  said  about  the  role  of  the  ionosphere 
in  VLF  and  ELF  radio  wave  propagation.  It  is  apparently  permis¬ 
sible  to  ignore  the  presence  of  the  ionosphere  at  distances  less 
than  about  30-50  kilometers,  in  considering  the  fields  near  the 
ground.^-  Extension  of  the  results  of  this  paper  to  such  dis¬ 
tances,  by  numerical  integration  of  the  branch  line  integral, 

would  therefore  enable  comparison  of  these  results  with  those 

4  _  K 

obtained  by  R.  K.  Moore  and  by  Banos  and  Wesley.  Their 
analyses  were  concerned  with  submerged  antennas  when  the  atmos¬ 
phere  is  assumed  homogeneous . 

4 

R.  K.  Moore,  Doctoral  Dissertation,  The  Theory  of  Radio 
Communications  Between  Submerged  Submarines,  Cornell  llniv.  1951 . 

5  ~ 

A .  Banos,  Jr.,  and  j.  P.  Wesley,  The  Horizontal  Dipole  in 
a  Conducting  Half-Space,  University  of  California  Marine 
Physical  Laboratory,  Reports  No.  53-53,  1953,  and  No.  5^-33,  195^. 


2.0  THE  SUBMERGED  VERTICAL  ELECTRIC  DIPOLE 

2.1  Introduction  of  Equations 

The  development  is  in  terms  of  the  Hertz  vector,  n,  from 
which  the  electric  and  magnetic  fields  are  derivable  by  the  re¬ 
lations: 

4 

f  =  M.sv  x  (2.1-1) 

f  =  vv.n  -  iii  2-H.  (2.1-2) 

atz 

where  e  =  ( e  +  ^j)  is  the  so-called  complex  dielectric  constant, 

-4 

and  I!  satisfies  the  vector  equation 

2"> 

v  X  v  X  n  -  vv-n  +  lie  ^-4  =  0  (2.1-3) 

dt 

The  constant  e  is  given  by  keQ,  where  k  is  the  ordinary  dielec- 
trie  constant,  and  eQ  is  given  by  IO'/^ttc  ,  with  c  =  velocity  of 
light  in  free  space.  Conductivity  is  represented  by  a,  while 
od  is  2irf. 

These  equations  are  applicable  in  regions  which,  indi¬ 
vidually  are  homogeneous  and  isotropic,  and  in  which  no  free 
charges  are  present.  In  the  remainder  of  our  discussion  it  will 
be  assumed  that  n(x,y,z,t)  =  n(x,y,z)  e'^ut.  The  time  factor 
eitDt  will  be  suppressed,  and  the  Hertz  vector  will  be  understood 
to  be  a  function  of  position  only.  With  this  modification, 
equation  (2.1-3)  can  be  written: 

vxvxn-vv-n  -  k2n  =  0  ( 2.1-4) 


-7- 


-8- 


2—22  2 
where  k  =  (ifico  =  txeco  -  ieup,a .  If  the  operator  v  is  under¬ 
stood  as  acting  on  rectangular  components  of  if,  we  use  a  common 
vector  identity  to  simplify  further: 

v2  if  +  k2  Tf=  0  (2.1-5) 

Thus  the  rectangular  components  of  if  individually  satisfy  the 
homogeneous  equation,  in  source-free  regions.  To  take  account  of 
the  source,  it  is  necessary  only  to  add  solutions  of  (2.1-5)  to 
a  solution  of  the  inhomogeneous  wave  equation,  in  such  a  way 
that  the  vector  is  singular  at  the  source  in  the  prescribed 
manner.  In  the  present  problem  this  is  accomplished  by  designa¬ 
ting  (see  Figure  2.1)  for  medium  1: 

-  S(pri)  +  ?1  (sac)  f2-1'6' 

where  satisfies  the  inhomogeneous  wave  equation,  and 

if,  sec)  the  homogeneous  wave  equation,  ifg  and  lf^  have  only  to 
satisfy  the  homogeneous  equations  as  written  for  their  respective 
media: 

n  =  2,  5  (2.1-7) 


v2  If  +  k2  If  =  0 
n  n  n 


In  order  that  represent  properly  the  field  of  the  electric 

dipole  source  oriented  along  the  z-axis,  it  is  written  specifi¬ 
cally  as  : 

-  _-ik,R,  . 

it  (2.1-8) 


#  I  dl  e"iklRl 

n(Pri)  ~  - KT“ 


4ir  inis'  *'1 

with  ic  the  unit  vector  in  the  z-direction.  In  cylindrical 

I  o  2 

coordinates  obviously  =  V p  +(z+zQ)  .  The  problem  has  cylin¬ 
drical  symmetry,  and  is  subject  to  the  boundary  conditions: 
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Figure  2 . 1 
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z  =  0  ic  x  (^2  -  £]_)  =  0 

it  X  (i!2  -  3X)  =  0  (2.1-9) 

z  =  h  it  x  -  ?2)  =  0 

it  x  -  32)  =  o 

where  the  numerical  subscripts  1,2,3  denote  fields  in  media 
1,2,3,  respectively.  An  additional  boundary  condition  is  the 
so-called  "radiation  condition" s 


lim 
R  -*■  00 


R  (gj  +  ik*)  =  0 


(2.1-10) 


where  R  is  radial  distance  and  f  is  any  of  the  field  components. 

We  now  assume  that  the  fields  everywhere  can  be  represented 
in  terms  of  the  z-axis  oriented  Hertzian  potential;  if  the  final 
solution  satisfies  the  boundary  conditions  and  behaves  properly 
as  -*•  0,  the  assumption  will  have  been  justified.  Consequently, 
we  henceforth  drop  the  vector  designation,  and  write  simply 
"II",  with  the  orientation  understood.  In  cylindrical  coordinates 
we  have; 


1  s 
p  37 


(p 


an 


3p“ 


£) 


i  a2n„  o 

+  K  — 7—  +  k  2n 
p7  30  n  n 


=  o 


(2.1-11) 


One  expects  that  solutions  may  be  represented  as  superpositions 
(discrete  or  continuous)  of  cylindrical  wave  functions; 


tn(P,0,z) 


Up) 

Up) 


(2.1-12) 


where  \  and  m  are  separation  parameters . 
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Since  the  solutions  are  required  to  possess  azimuthal 
symmetry,  and  must  be  finite  at  p  =  0,  we  choose  the  Bessel 
functions  of  the  first  kind  and  m=0.  Thus  the  desired  func¬ 
tions  ares 

JQ(Xp)  e*  ^An  “  ^  (2.1-1J) 

2  2  2 

It  will  be  convenient  to  use  8  ,  defined  by  8  =  k  -  X  ,  as 

n  ■*  n  n  ' 

a  new  separation  parameter. 

In  the  three  regions  we  may  now  anticipate  the  representa¬ 
tions: 

—  R  o 

ni  =  M  {e  +Jb1(x)J0(Xp)e±  i^lz  dx}  (2.1-14) 

n2  =  M  J*  [a2( X)  e-i^2z  +  b2(x)  ei^2z]  JQ(xp)dX  (2.1-15) 
n5=Mja5(x)  e_i^3Z  J0(Ap)dX  (2.1-16) 

where  the  spectral  densities  b^,  b2,  a^,  a2  are  to  be  determined 
from  the  boundary  conditions,  and  the  limits  of  integration  will 
be  specified  later.  M  is  simply  an  amplitude  factor  depending 
on  the  "strength"  of  the  antenna,  and  is  described  by 

M  =  1  61~  (2.1-17) 

47ria)e 

where  I  dl  represents  the  dipole  moment  of  the  antenna. 

In  the  integrals  of  (2.1-14)  and  (2,1-16)  a  boundary  condi¬ 
tion  has  already  been  applied;  namely,  that  for  z  <  -zQ  and 
z  >  h,  only  "outgoing"  waves  are  physically  allowable.  We  note 
in  this  connection  that  the  sign  of  the  imaginary  part  of  jS  is 
to  be  taken  as  negative,  in  order  that  the  integrals  converge. 


i 
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2  2  2 

We  have  already  defined  k  =  <x>  p,e  -  io^ia,  in  which  ao  |ie 

and  cujxa  are  positive  real  numbers.  In  computing  k  we  will  take 
the  square  root  which  has  a  positive  real  (therefore  negative 
imaginary)  part.  It  can  then  be  shewn  that  the  primary  excitation 
can  be  represented  as  follows s 


e-lklRl 


OO 


Jo^  e  1 


-i£.  \z+zr 


dX 


(2.1-18) 


It  may  be  noted  from  figure  2.1  that  in  the  exponent  under  the 
integral,  zq  is  a  positive  number.  From  (2.1-18)  we  will  surmise 
that  the  limits  of  integration  in  (2.1-14),  (2.1-15)  and  (2.1-16) 
are  also  to  be  taken  from  zero  to  infinity.  For  region  1  we  have, 
finally: 

n1  =  M  Io  r  e-i^l  I  z+zo  I  +  bx(x)  e^l2  ]  Jo(Xp)dx 

(2.1-19) 

We  now  apply  the  boundary  conditions.  The  tangential  com¬ 
ponents  of  it  and  it  are  required  to  be  continuous  on  the  planes 
z  =  0  and  z  =  h.  it  and  3  are  derived  by  the  equations 


*n  =  VV*  *n  +  n  ?n 

it  =  T” -  7X 

n  icD(j,n  n 


(2.1-20) 

(2.1-21) 


and  are  found  to  be 


a2n 


E 


n 


-  •  E  n  =  0  *>  E 

np  3p8z  n?  nz 


(f- 


+  k. 


D  “n 


H, 


np 


H 


kn  9nn 

n  n 


SB  '  P 

h  0  ici4in  3  P  }  nz 


( 2.1-22) 


0 


-13- 


An  integration  of  E  and  can  be  performed  from  p  to 
infinity.  Since  the  fields  and  their  derivatives  are  required 
to  vanish  at  infinity,  the  boundary  conditions  are  then  speci¬ 
fied  in  a  somewhat  simpler  manner  by 


ann+l 

n=l, 

z=0 

m 

9  z 

9  z 

n=2 , 

z=h 

k2n  = 

n  n 

kn+l  nn+l 

n=l, 
n=2 , 

OXi 

II  II 
N  N 

in  which  we  have  assumed  nn  =  (iQ,  the  magnetic  permeability  of 

free  space,  for  all  three  media. 

The  conditions  (2.1-23)  are  applied  to  equations  (2.1-19)> 

(2.1-15)  and  (2.1-16).  The  algebra  is  somewhat  simplified  if 

the  coefficients  b^,  a2,  b2  an<5  are  assumed  to  contain  a 

factor  — - —  e^lZo  ,  thus: 

^1 


b^(A)  -  B1(A)  e  ^lZo 

a2(A)  =  A2(A)  jjg—  e  i/3lZo 

(2.1-24) 

b2(A)  =  B2(A)  e“i/3lzo 

a^(A)  =  A^(A)  jjj—  e  1/3  lZo 

Four  equations  are  obtained  in  the  four  unknowns: 

MV1)  =  MB2^2> 

P2(B2  e^2h  -  -  -^3 

^1  (®i+l)  =  ^2  (®2"*"^2  ^ 


(2.1-25) 
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k^(B2  e^2h+  A,  e_1^2h)  =  k*  A,  e’^3h 


3  "3 


(2.1-25) 


If  we  define  quantities  Kn(X)  such  that 
ij3 

Kn(M  "  TT 


(2.1-26) 


where  P „  =  cr  +  ia>e„,  we  can  rewrite  these  equations  as 
n  n  n 

P1K1  (Bl-1)  ■  P2K2 (B2  "a2  ) 

P2K2  (B2  e^2h  -  A2  e-1^2h )  =  -  Pycy^  e'1^11 
P1(B1+1)  =  P2(B2+A2) 

P2(B2  e^2h  +  A2  e-i^2h )  =  PyVyT1^ 

It  is  then  readily  determined  that 


(2.1-27) 


-  K2(K1~  K3>  +  1  (K1K3~  k2 )  tan  ^2h 
1  K2  (k1+  K^)  +  i  (K1K5+  k|)  tan  £2h 


(2.1-28) 


Then  Kn(>)  can  be  considered  as  intrinsic  wave  impedances  for 
the  media  concerned.  It  should  be  pointed  out  that  rewriting 


B^  in  the  form 


+  K^tanh  ijSgh)  -  K2  (K?+  I^tanh  ip>2h) 

1  Kx(K2+  tanh  i£2h)  +  K^Ky-  K2  tanh  i^h) 


leads  quickly  to  the  representation 
K,  -  Z~ 


B1  = 


K1  +  Z2 


where  Z2  =  K2 


Kj  +  K2  tanh  i02h 
K2  +  tanh  i/32h 


(2.1-29) 

(2.1-30) 

(2.1-31) 
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Z2  is  the  wave  impedance  looking  into  medium  2  at  the 
boundary  between  1  and  2 . 

The  other  coefficients  are  found  to  be 

Pi  k7  +  z7 

A2^  =  ^  K~+~Z^ 

P-,  K9  "  Z9 

B2^^  K-rz^  (2.1-22) 

A-j (>> )  =  2K,  ^1  (K2  +  Z2^e_1^2h  ei^2h 

3  3  iSJ  (Kx  +  Z2)  (kJ  +  K5)  * 

The  integrals  for  the  Hertz  vector  in  each  region  now  become: 


II1  =  M 


I  z+2o  I  +  ^2  (2.1-53) 

o  1+  2  Ip7~ 


n2  =  2m 


f^l  K2cosh(i^2z)-Z2sinh(i^z)  ^  ^  ^ 

TW~  k  ■+•  7._  e  1  o  J (Ap)AdA 


o1^2 


K1  +  Z2 


n  _ _ 2m_ 

5  a,+ia>e- 


p®  (k2+z2) 
(kl+z2)  (k2+k5) 


(2.1-24) 

!_1^lZo  e_1^2h  e_1^2^2_h^J0(^p  )AdX 

(2.1-25) 


2.2  The  Hertz  Potential  in  Region  1: 
We  will  first  consider  region  1: 


o® 


n1  =  m 


G(jBr  B2,  03)  J0(Ap)AdX  (2.2-1) 

o  ^ 

where  G  (j3^,  jS2»  jS^)  is  1/ip^  times  the  bracketed  part  of  the 
integrand  in  equation  (2.1-22).  We  introduce  the  identity: 


J0(Ap)  =  1/2  Lh^  (Ap)  +  H2q  (Xp)J  (2.2-2) 

and  designate 

nx  =  (i1  +  i2) 


(2.2-2) 
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where 

co 

X1  “  I  JHo  (Xp)  G(^1'  ^2 '  Pj)  XdX  (2 .2 -4a) 

c0 

I2  =  I  J Ho  (xp)  G(/3i'  /V  ^  XdX  (2 .2 -4b) 

The  functions  in  this  and  the  other  integrands  occurring 
in  our  problem  are  analytic  except  at  isolated  singularities 
(poles  and  branch  points),  and  can  be  continued  into  the  complex 
plane  when  proper  account  is  taken  of  such  singularities. 

We  therefore  transform  our  integrals  along  the  real  axis 
of  into  contour  integrals  in  the  right  half  of  the  complex 
A-plane.  In  order  that  this  be  feasible,  it  is  necessary  that 
the  integrands  behave  properly  when  the  modulus  of  X  becomes 
large.  We  note  that  as  |  A  |  -CP  , 

6  ~  i  X  . 

n 

Also,  K  (x)~  A  X,  n=l,2,3  ,  where  the  A  are  constants.  The 
n '  n  n 

main  question  now  concerns  the  hyperbolic  tangent,  which  has 
poles  where 

cosh  (i^h)  =  0. 

These  occur  for 

i/32h  =  t(n+  j)  iir  ,  n  =  0,  1,  2,  5  -----  , 
which  requires 

X  =  t  J  k^  -  (n+  £)2,  n  =  0,  1,  2,  5  . 

2 

Since  k2  has  a  negative  imaginary  part,  the  poles  of  concern 
will  be  to  the  right  of  the  negative  imaginary  axis. 


o 
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The  customary  technique  for  dealing  with  these  poles  is  to 
restrict  the  contours  to  a  set  of  circular  arcs  or  rectangles 
such  that  the  contours  always  pass  between  the  poles.  The  hyper¬ 
bolic  tangent  is  analytic  everywhere  except  at  such  poles,  so  it 
is  bounded  on  this  set  of  contours . 

Therefore  as  |  X  |  -*oo  ,  on  the  appropriate  set  of  contours, 

Z2  ~  A2^ 

where  A2  is  a  constant.  Then 


K 

K 


2 

2 


(2.2-5) 


The  asymptotic  behavior  of  the  Hankel  functions  is : 

M  *  i(x’,”r/4)  (2-2^> 

„2  M  _  fli ~  e-ifV-  rA)  (2. 2 -6b) 

of  which  the  first  is  valid  when  -tt<  arg(Xp)<  2tt,  and  the  second 
when  -  2tt  <  arg(Xp)  <  tt. 

This  behavior  of  the  Hankel  functions  ensures  that  if  the 

imaginary  part  of  \  is  positive,  the  integrand  of  1^  vanishes 

faster  than  1  as  |X|  -*  o»  .  If  the  imaginary  part  of  X  is 

1 

T~\as 

|X  |  -*  eo  ,  on  the  appropriate  contours.  We  will  therefore 
transform  1^  into  an  integral  in  the  first  quadrant  and  ^  into 
an  integral  in  the  fourth  quadrant. 

It  is  safe  to  assume,  without  going  into  delicacies  of 
mathematical  proof,  that  no  poles  of  the  integrand  can  occur 
in  the  first  quadrant.  If  such  poles  (X  )  did  occur,  they  would 


negative,  the  integrand  of  ^  vanishes  faster  than 
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give  rise  to  modes  of  propagation  for  which,  at  large  distances, 
the  field  expressions  would  contain  exponentials  of  the  type 
elAhP .  Since  the  time  factor  is  ela>t,  such  expressions  would 
represent  fields  traveling  in  the  negative  p  direction,  that  is, 
waves  converging  onto  the  source.  The  radiation  condition  pre¬ 
cludes  the  possibility  of  such  waves. 

We  use  this  assumption  and  the  estimates  of  the  integral 
1^  as  |>»  |  -oo  to  deform  1^  into  the  path 


T  M 

*1  “  7  J 


h„  (Xp)  <3(0!,  02,  03)  xax  . 


Using  a  real  parameter  of  integration  u,  with  A 


=  1U, 


i  -  M 

Ii  ~  '  7 


Ho  (iup)  6(/3lf  A2,  A5)  udu 


=  (“)  ^  J  Ko(up)  6(jB1,  /32,  /3?)  udu 


(2.2-7) 

(2 .2 -8a) 
(2 .2  -8b) 


where  KQ(up)  is  defined  in  the  conventional  way: 


Kq (x )  =  (ix)  . 

In  the  fourth  quadrant  we  expect  to  find  poles  and  branch 
points.  The  branch  points  arise  from  the  quantities  ^n=V  . 

Let  us  now  treat  the  I2  portion  of  (2.2-3): 


X2  =  7  j  Ho  G^l'  ^2 '  ^  AdA  ' 

In  Figure  2.2  is  a  representation  of  the  complex  X-plane, 
on  which  branch  cuts  have  been  drawn.  It  is  convenient  in  this 
case  to  let  them  be  vertical  downward;  if  any  of  these  branch 
lines  happens  to  cross  a  pole,  it  is  considered  to  be  deformed 
slightly  to  avoid  the  pole. 
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Figure  2.2.  Path  of  Integration  in  4th  Quadrant. 


I 
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We  have  already  verified  that  integration  along  the  appro¬ 
priate  infinite  quarter  circle,  or  portions  thereof,  contributes 
nothing.  The  determination  of  I2  is  thereby  reduced  to  an  inte¬ 
gration  along  the  negative  imaginary  axis,  integrations  looping 
back  around  each  branch  cut,  and  an  evaluation  of  the  residues 
at  any  poles.  Thus, 

_  f  M  f  'J\*>  M  P  M  f  M  P  /Mx  y  .  V  Roa  ? 

I2  "  [  7  Jo  +  7  J  +7  J  +  7  J  ~  (7)  2iri  L  Res*  j 

B.L.l  B.L.2  B.L.3 

(2.2-9) 

The  first  integral  in  the  brackets  of  (2.2-9)  is ' 

-i  CD  © 

“  J  H2(Ap)G(j3lf|32,/33)Xd^  =  -  |  jH2(-iu0)6(/B1,/32,|35)udu 

CO 

=  “(7)^7  J  K£>(up)G(|B1,i32f/35)udu 

(2.2-10) 

This  will  exactly  cancel  1^,  as  given  by  (2. 2 -8b). 

We  note  now  by  reference  to  (2.1-55)  and  the  definition  of 
Z2  that  the  integrand  is  single-valued  on  the  branch  cut  for  k2« 

The  integration  upward  along  the  left  side  of  this  branch  cut 
consequently  cancels  the  integration  downward  along  the  other 
side,  and  the  problem,  in  its  entirety,  is  further  reduced  to 
integration  along  branch  cuts  2  and  5,  and  the  sum  of  residues. 

We  will  now  list  the  formal  expressions  for  the  integrals 
along  branch  lines  2  and  5*  The  discussion  of  evaluating  these 
is  postponed  until  the  next  section,  in  which  the  magnitudes  of 
the  various  quantities  involved  are  taken  into  account  and  suit¬ 
able  approximations  made.  Without  such  approximations  it  is 

1 4 

,  t 
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more  or  less  impossible  to  deal  with  the  integrals  in  a  general 
way. 

k. 


r  M 

CBL2  =  7 


p? 


k,-i  co 

r  , 


Ho(kp)G01^2^5)^dX  +  |  J  H‘(Xp)G01,A2,A5)XdX 


k^-i  Oo 


M 

7 


k, 

p3 


h2(>.p)[g(?1,p2,p  )-a(^,p2,p  )]  XdX  (2.2-11) 


kj-i  » 


is  J~ 


2  2 

where  (3^  is  V  k^  -  X  evaluated  as  one  approaches  the  branch 
cut  from  the  left.  Likewise, 
k 


M 


BL3  -  7  J 


1 

1  H^(XP )  )-G(-^1,/32,/3  )]  XdX  (2.2-12) 


k^-i  co 


/  9  2 

where  ^3 is  v  k^-X  evaluated  as  one  approaches  this  branch  cut 
from  the  left. 

It  remains  to  specify  the  contribution  from  residues  at 
poles.  Since  the  evaluation  is  at  (K^+Zj )  =  0,  the  integrand 


concerned  becomes 


0  K.  -  Z0  ifi,  (z-z  ) 

H^(Xp)  ,  J-  2 - It - —  *  . 


K1  +  Z2 


(2.2-15) 


If  the  poles  are  designated  by  Xn#  the  formal  expression 
for  the  contribution  from  the  residues  is  now  (assuming  the 
poles  are  simple): 

-2*1  («)  £  Res  >  -  l  H2(V)  (2.2-14) 

n 
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where  fl?  =  V  k?  -  X^  .  The  X  are  the  zeroes  of  K,  +  Z0.  The 
rl  Inn  1  1 

matter  of  finding  these  zeroes  is  deferred  until  the  next  section, 
where  specialization  to  the  water -air-ionosphere  problem  is  made. 


2. 3  The  Hertz  Potential  in  Region  2. 

The  same  procedure  is  followed  for  the  development  in  this 
case.  The  integral  concerned  is  (36).  In  examining  the  behavior 
of  this  integrand  for  large  |X  |  ,  we  note  that  the  bracketed  term 
approaches 


cos  (-iXz  ) -1A2X  sin(-iXz)  A£  e“ 

A^X  +~ A^"X  A^  +  A2 

where  A^  and  are  constants.  For  ReCX]  >  0  and  z  &  0  this  term 
is  evidently  bounded  along  an  infinite  radius  in  the  right  half 
plane,  and  the  exponential  behavior  of  the  Hankel  functions  en¬ 
sures  that  the  integrands  vanish  faster  than  1/|X|  as|A|  -»  *  , 
as  before,  along  the  appropriate  contours. 

As  before  we  will  be  left  with  a  representation  involving 
integration  along  the  branch  cuts  and  residues  at  any  poles. 

Again  the  integrand  is  single  valued  around  the  branch  cut 
for  ^2  so  that  the  branch  line  1  integration  gives  zero.  The 
branch  line  2  integration  can  be  written: 


XBL2  =  _iM 


k, 

p  3  2 


H^(Xp)  Cf(/31,jS2,P3)  -F01,j32nS,)]  XdX  (2.3-1) 


k,-ico 

where 

-ifl,z_  IC,  cos  S0z-iZ0  sin  B0z 

rp  f  q  q  a  N  _  ® _ 1  O  t _ 2 _  ■£* 

F\Pi'P2'Py  ^  k1  +  z2 


(2.3-2) 
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hv  = 


The  branch  line  3  integration  is  represented  by 

ki 

hJ(X5)  Cf  01,?2,P3)  xax  (2.5-3) 


k^-i  oo 

The  contribution  from  residues  at  poles  is  given  by 

r  K2  cos  ^z-i  K]_  sin  i^2z 


_  o  -ifi^z  f  ] 

2ttM  T  H*(X  p)X  - — '  • 
Lt  o '  nH  7  n  an 

Po 


(Kx  +  Z2) 


1 


(2.3-4) 


x=x_ 


where  obviously  the  Xn  are  the  same  as  for  the  corresponding 
expression  in  region  1. 


The  Hertz  Potential  in  Region  3 

Since  for  our  problem  this  is  of  relatively  less  importance 
than  the  regions  1  and  2  potentials,  we  will  hereafter  omit  it 
from  our  discussion. 


2.4  Specialization  to  the  Water -Air-Ionosphere  Problem 

? 

In  the  1-1000  cps  range  of  frequency  k^  is  very  accurately 
given  by  the  term  (-iu*ia^).  This  is  equivalent  to  saying  that 
the  displacement  current  in  the  wave  equation  is  negligibly 
small  compared  with  the  conduction  current.  The  criterion  that 
this  be  true  is  that  c/W  be  much  greater  than  1. 

The  situation  is  not  so  clear-cut  in  the  case  of  the  iono¬ 
spheric  propagation  constant.  For  one  thing,  there  is  con¬ 
siderable  doubt  as  to  the  proper  value  of  conductivity  to  be 
used.  In  this  paper  we  are  assuming  a  physical  model  which 
represents  an  extremely  simplified  picture  of  the  actual  situa¬ 
tion.  The  ionosphere  is  not  homogeneous  and  not  sharply  bounded; 
moreover,  the  bottom  portion  of  its  lowest  ionized  region  has 


1 


1 
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diurnal  and  seasonal  variations.  Under  such  circumstances,  the 
most  that  can  be  expected  of  our  theoretical  development  is  that 
it  gives  answers  which  are  in  reasonably  good  correspondence  with 
observed  phenomena.  The  "correct"  values  to  be  used  for  a  and  h 
are  those  which  give  results  conforming  most  closely  with  obser¬ 
vation.  If  the  discrepancies  are  too  great,  then  another,  per¬ 
haps  more  elaborate  model  must  be  used. 

The  effective  conductivity  (the  conductivity  which,  for 
the  model  of  a  sharply  bounded  homogeneous  ionosphere,  gives 
results  that  compare  favorably  with  experiment),  can  be  estimated 
from  the  empirically  developed  Austin-Cohen  formula^  which  gives 
the  attenuation  in  field  strength  from  one  point  to  another  when 
transmitter  and  receiver  are  located  on  the  earth’s  surface. 

Such  a  comparison  is  made  by  Bremmer^,  who  then  concludes  that 

i  -4 

a  conductivity  of  around  4  x  10  mhos  per  meter  is  appropriate. 
-4  8 

The  value  of  10  is  also  used  by  Schumann.  Recent  studies 

seem  to  show,  however,  that  this  value  is  too  great.  A  large 

number  of  "atmospherics",  the  electromagnetic  waves  radiated  by 

lightning  discharges,  were  studied  over  a  period  of  several 

q  in 

years  by  Chapman  and  Edwards,-^  Chapman  and  Matthews, 

^Austin  and  Cohen,  Bulletin  of  the  Bureau  of  Standards, 

7,  315  (1911) 

7 

H.  Bremmer,  Terrestrial  Radio  Waves  (Elsevier  Publish¬ 
ing  Company,  Amsterdam  1949)  p.  22b 

®W.  0.  Schumann,  "Uber  die  oberfelder  bei  der  Ausbreitung 
langer,  elektrischer  We  lien"  Zeitschrift  fur  angewandte  Physik 

6,  35  (1959) 

q 

7F.W.  Chapman  and  A.G.  Edwards,  Proc.  U.R.S.I.  8,  pt.2, 

351.  (1950)  -  “ 

10F.W.  Chapman  and  W.D.  Matthews,  Nature  172,  495  (1953) 
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11  12 
Chapman  and  Jolley  and  Chapman  and  Macano.  Spectral 

analyses  were  performed  over  the  range  40  cps  to  16  kc/sec, 

and  marked  differences  were  found,  below  10  kc/sec,  from  the 

attenuation  to  be  expected  on  the  basis  of  the  Austin-Cohen 

13 

formula.  The  whole  matter  is  discussed  by  Wait,  ^  who  finds 
that  good  agreement  is  obtained  between  theoretical  and  experi¬ 
mental  results  for  VLF  if  the  ionospheric  conductivity  is  taken 

to  be  about  one  two-hundredth  of  the  value  used  by  Bremmer, 

-6 

i.e.,  i*i  2  x  10  mhos/meter.  In  the  ELF  range  there  are  in¬ 
dications  that  penetration  of  electromagnetic  energy  is  deeper 
into  the  ionosphere,  so  that  the  effective  conductivity  as  well 
as  the  effective  height  of  reflection  is  greater  than  the  in¬ 
dicated  value  for  VLF.  A  recent  paper  by  Wait1  discusses  this 
matter. 

In  the  frequency  range  100  cycles/sec  to  1000  cycles/sec 

there  is  some  indication  that  the  attenuation  rates  are  given 

more  exactly  when  a  more  elaborate  model  of  the  ionosphere  is 

used  in  the  theoretical  development.  An  analysis  is  made  by 

14 

Wait  for  an  atmosphere -ionosphere  model  wherein  the  refractive 
index  is  given  as  a  function  of  height  by: 


1  F.  W.  Chapman  and  A.  Jolley,  Proc.  U.R.S.I.  10,  part  2, 
355  (1950) 

12 

F.  W.  Chapman  and  R.C.V.  Macano,  Nature  177/  931/ 

(1956) 

^J.  R.  Wait  and  A.  Murphy,  Proc.  I.R.E.  45  ,  754  (1957 
^J.  R.  Wait,  Journal  of  Geophys.  Research  65/  597  (I960) 
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N(z)  = 
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l 

N  ez_h//L 


0  <  z  <  h 
h  <  z 


(2.4-1) 


In  these  equations  h  is  the  height  at  which  the  refractive  index 

changes  discontinuous ly  from  1.0  to  some  value  N,  and  L  is  a 

scale  factor.  Best  agreement  with  the  results  of  Chapman  and 
12 

Macario  seems  to  be  obtained  when  L  is  chosen  as  being  equal 

to  JO.  The  result  is  that  the  attenuation  in  nepers  per  unit 

1/2 

length  does  not  vary  as  f  ,  as  it  would  for  a  homogeneous 

atmosphere,  but  as  f\  1/2  <  y  <  1. 

The  experimental  data  in  the  ELF  range  is  open  to  a  certain 

amount  of  question,  however,  and  for  the  submarine  problem  it 

would  seem  most  logical  to  proceed  with  an  analysis  based  on  the 

homogeneous  ionosphere.  If  necessary,  a  slight  modification  may 

possibly  be  made  in  the  final  answer  to  account  for,  at  least, 

the  gross  effects  of  an  inhomogeneous  ionosphere.  We  will  take 

4 

the  height  of  the  lower  boundary  of  the  ionosphere  as  9  x  10 
meters,  and  a  conductivity  of  10-^  mhos/meter,  these  figures 
corresponding  to  ones  which  have  been  used  by  Wait  in  ELF 
studies . ^ 


To  go  back,  then,  to  the  criterion  o/coe  »  1,  we  see  that 
the  conduction  current  is  much  greater  than  the  displacement 
current  provided  f  «  1.8  x  10^  c.  This  is  quite  well  satis¬ 
fied  for  the  ELF  range,  and  will  be  assumed  throughout  the  re¬ 
mainder  of  our  discussion. 
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In  air,  the  conduction  current  is  negligibly  small  compared 

ic 

with  the  displacement  current.  J  For  the  three  regions,  then: 


k*  * 

ic^oCTl 

k*  * 

A0eo 

(2.4-2) 

k2 

Kj  S3 

-  ionoa5  . 

Therefore 

KX(A) 

i/3  ^ 

k2(a) 

ft  2 

**  (DC 

(2.4-5) 

k^(a) 

O 

* 

2.5  The  Hertz 

Potential  in  Water 

The  branch  line  2  integral  for  region  1  will 

now  be  in- 

vestigated  for 

the  water -air -ionosphere  problem. 

From  the  re- 

suits  of  the  previous  section  and  using  (2.1-28)  or  (2.1-29) 
one  gets 

til (q.2  — ol jT )  —  a-2  (a-j+c^T) 

B1  "  a1(a2-a^T)  +  (^(aj+agT) 


a^c^-a^T)  -  a^(a2+a-LT) 
a^c^+c^T)  +  a^c^-a^T) 


(2.5-1) 


15  2 

-Although  k2  is  practically  real,  it  is  nevertheless  im¬ 
portant  in  arguing  the  validity  of  the  original  formulation  of 
our  solutions  that  it  always  be  understood  to  retain  a  small 
negative  imaginary  part;  hence  it  is  not  on  the  real  axis  of  A. 
This  matter  is  discussed  thoroughly  by  G.  M.  Wing,  Mathematical 
Aspects  o £  the  Problem  of  Acoustic  Waves  in  a  Plane  Stratified 
Medium.  Sandia  Corporation  Memorandum,  i960. 
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where 


_  ^1 
1  "  °1 


a2  =  £7^  '  a3  =  '  T  =  tanA2h'  ( 2  •  5 -2 ) 


and  B^  has  its  meaning  as  previously  specified.  By  reference 
to  equations  (2.1-53)  and  (2.2-1)  we  see  that 

i/3.(z -z  )  _ 

g  » ^2 ' ^3 )  (£ ]_ » ^2 '  ”^3 ^  =  ip  Lb ^ (a^» Q2 * ® ^ ) (a » ^2 ' ^3 ^ 


The  bracketed  part  of  this  can  be  written 


where 


Q-ouP  Q+a,P 

_ 12 _ 2_ 

R+a-,S  R-a,S 

3  3 


Q  =  (ai~a2T) 
P  =  o^+djT 

R  =  a2  (ctj+c^T) 


2a^  (SQ+PR) 

"T2 - T~7 

R  -  a-,  S 

3 


(2.5-3) 


S  =  a2  -a^T 


This  becomes 


-  4a1a2^a^ 


— 2 - 2 - T - 7 — 

a2  (djC+c^s)  —  d^  (d2^— ci^s) 

where  S  =  sin  ^h,  c  =  cos  It  will  be  convenient  to 

rewrite  this  as 


D(X) 

a. 


(2.5-4) 


where  6.  =  — 

1  CD€, 


6,  =  — ,  and 

3  0)€o 


D(X)  =  b^22{^1C+6^2S)2  -  -  6^2C)2 
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The  branch  line  2  integration  is  now 


hoi  -  2U44165J 


e  l'Sl^z_zO^H2  (Xpjxax  (2.5-5) 

k, -i»o 
p 

We  will  designate  k^  =  a(l-i)  and  change  the  parameter  of  integra¬ 
tion  by  writing  X  =  a(l  -i  -it).  It  will  be  necessary  to  use 
the  asymptotic  approximation  to  the  Hankel  function,  (42),  which 
is  quite  good  for  large  values  of  the  argument  and  which,  at 
least  when  the  index  is  zero,  is  within  an  order  of  magnitude  of 
the  correct  value  down  to  values  of  the  argument  only  a  little 
larger  than  one.  (See,  for  example,  Jahnke  and  Emde,  Tables  of 
Functions,  Dover,  1945)*  In  terms  of  distance  and  frequency 
this  requires,  roughly,  p'/f  &  1.2  x  10^. 

The  branch  line  integral  is  then 


Irt9  =  -2Ma5/2  ei7r/4  6.6,  JH  ft1  FOOe^dt 

BLZ  1  P  TTp  'fc 


(2.5-6) 


where 

■J t+2  (l+i)  (1-i-it)1^2  g1/*!  (z  -z0) 

F(t)  * - Errrtn - 

in  which  the  j3's  are  expressed  as  functions  of  t.  F(t)  is 
analytic  at  t  =  0  and  is  bounded  on  the  real  axis  of  t.  It 
follows  from  Watson's  lemma1^  that  an  asymptotic  expansion 
exists  and  can  be  obtained  by  integrating  the  terms  resulting 
from  a  Taylor  series  expansion  of  F(t)  around  t  =  0. 


^Jeffreys  and  Jeffreys, 
(Cambridge  University  Press, 


Methods  of  Mathematical  Physics, 
1946)  pp.  471,  668. 


-jo¬ 


lt  is  readily  found  that 

2  eibl<2-zo) 


F<0)  »  — 7 


[t^cos  hb2+  6^2  sin  hbj] 


(2.5-7) 


where 


bx  =  V  kj2  +  2ia2 

b2  =  ^2^  +  2  ia^  . 

2  2 

Since  k^  »  2ia  ,  it  will  be  sufficiently  accurate  to  take 

bi  "  kr 

The  first  term  arising  from  an  integration  over  the  Taylor 
series  expansion  of  F(t)  is  therefore  proportional  to: 


OO 

F(°)  v't  e"aptdt  =  j  F(0)|-a^- 


J/2 


(2.5-8) 


If  the  second  term  were  obtained,  it  would  contain  the  factor 


F'  (0) 


(2.5-9) 


and  so  forth,  in  increasing  inverse  powers  of  p.  With  multi- 

-1/2 

plication  by  the  factor  p  from  outside  the  integral,  it  can 

be  concluded  that  the  branch  line  integral  is  attenuated  at 

least  as  e~aD/p^  at  large  distances.  Since  the  factors  whose 

product  is  F(t)  are  smoothly  varying  near  t  =  0,  it  is  expected 

that  F'(0),  F"(0),  etc.  will  not  be  excessively  large  compared 

with  F(0),  and  that  as  an  estimate  when  p  is  of  the  order  of 

1000  km  or  greater,  it  will  suffice  to  consider  only  this 
-2 

p  term. 

The  complete  expression  for  the  branch  line  integral  is: 
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«  _  (z -z  ) 

T  o  «■  M  -iJir/4  6la  e  1 

It,t?  <u  i.*JL  M  e  - yr 

65(b10+61b2sr 


e-ap(l+i) 

- 2 - 

P 


(2.5-10) 


where  we  have  abbreviated  C  =  cos  hb2  and  S  =  sin  hb2 . 

When  |k2 |  «  |kj|,  which  requires 

f  «  2  x  1010  c^  (2.5-11) 

the  form  of  F(t)  can  be  simplified  somewhat.  Under  this  con¬ 
dition  D (t)  can  be  written 


D(t)  *  j924(t)S2(t)  .  (2.5-12) 

Then  (2.5-5)  assumes  the  form 

ei/3l(Z-Zo>  — Hq2  (Ap)AdA  (2.5-13) 

j  -  £o  s 

k_. -ico 

3 

Substituting  A  =  a(l-i-it)  and  using  the  asymptotic  expression 
for  the  Hankel  function, 


f  2iM 

JBL2  *  TZT. 


BL2 


-  -2*telr/l‘  -4?  Stir 


where 


F(t)  =  e^l<z-zo) 

^2  s 


Evidently 

F(o)  -  -fj. 
b2  5 

Making  the  expansion  of  F(t)  in  Taylor  series  therefore  gives 
as  the  first  term  estimate  of 

*  ..  2  J7e'iJnr/ka  _  ikn  (z  -z  )  e-a^1+i) 

xBL2  m  m  777-2=7 -  e  1  0  - 2 -  (2.5-15) 

6l63b2  s  p 
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The  integral  around  branch  cut  3  must  now  be  considered. 

One  can  see  rather  quickly,  however,  that  the  contribution  from 
this  integral  must  be  quite  small  at  large  distances.  This  is 
so  because  of  the  exponential  term  arising  from  the  Hankel 
function.  For  A  in  the  vicinity  of  k^,  which  has  an  imaginary 
part  equal  to  about  3-98  x  10""^  and  for  p  around  100  km, 
obviously  the  attenuation  will  be  enormous  for  ELF  or  higher 
frequencies.  We  will  therefore  neglect  this  contribution. 

The  next  task  is  to  solve  the  equation  +  Z^  =  0,  in 
order  that  evaluation  of  the  expression  (2.2-14)  can  be  made. 

The  roots  An  cannot  be  determined  generally,  and  it  is  necessary 
to  consider  our  situation  in  a  rather  specific  way. 

From  the  values  for  Kn(X)  given  by  (2.4-3 )•  it  can  be 
determined  that  the  zeroes  of  K1  +  Z^  are  obtainable  by  solving 
the  equation 


tan  h^2  ^  fij03  + 

~W  '  -  ^?22 


(2.5-16) 


This  is  quite  a  formidable  task.  It  is  useless  to  search 
for  real  values  of  A  ,  since  this  would  imply  lossless  propaga¬ 
tion  in  the  radial  direction.  In  the  analysis  of  a  physical 
problem  such  propagation  will  appear  to  exist  only  in  the  case 
when  some  non -realizable  perfection  has  been  assigned  to  the 
media  in  which  propagation  occurs,  i.e.,  perfectly  reflecting 
boundaries,  and/or  perfect  --  that  is,  lossless  --  dielectrics. 
Since  no  such  assumptions  have  been  made  in  this  problem,  we 
immediately  resign  ourselves  to  the  requirement  that  the  roots 
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>*n  will  be  complex.  This  can,  in  fact,  be  verified  from 
(2.5-16)  without  great  difficulty. 

At  this  point  it  is  possible  to  invoke  a  physical  argument. 

If  the  earth  and  ionosphere  boundaries  were  perfectly  conducting, 
only  the  so-called  TEM  mode  of  propagation  could  have  a  real 
propagation  constant,  for  the  1-1000  cps  frequency  range.  This 
would  be  exactly  ^2-  The  propagation  constants  for  other  modes, 
if  any  were  excited,  would  be  purely  imaginary  and  such  modes 
would  be  exponentially  damped  to  a  very  high  degree.  This  is 
because  the  distance  between  boundaries  is  considerably  less  than 
one  half  wave  length  for  any  frequency  in  this  range.  The  intro¬ 
duction  of  finite  conductivity  in  the  walls  causes  nonzero 
electric  field  tangential  to  the  walls  but  as  long  as  the  con¬ 
ductivity  can  be  considered  high  (a/ cue  »  1)  the  propagation  is 
still  considered  to  be  of  the  TEM  type.  As  the  resistivity 
varies  continuously  from  zero  to  a  small  but  non-zero  value, 
the  associated  propagation  constant  must  vary  continuously  from 
k2  to  some  value  in  the  vicinity  of  k2 .  If  any  poles  were  on 
the  imaginary  axis,  they  would  also  move  continuously  to  new 
positions,  but  their  new  locations  would  not  be  drastically 
closer  to  the  real  axis.  This  matter  is  discussed  in  detail  in 
Appendix  A.  On  the  basis  of  the  results  developed  there,  we 
venture  to  assert  that  only  one  root  will  be  found  in  the  vicinity 
of  the  real  axis,  and  that  this  root  will  not  be  too  far  removed 
from  k2 .  We  will  therefore  use  the  small  argument  approximation 
to  the  tangent: 


tan  h^  «  h^ 


(2.5-17) 
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which  is  reasonably  valid  if  IhjSjl  5  0.4,  implying 

T 


7k  „2 


-  V 


<  4.43  x  10' 


4 

where  an  ionosphere  height  of  9  x  10  meters  has  been  assigned. 

2  2 

Letting  the  root  be  ?^o,  we  represent  X*  as  k£  (1  +  A);  then 


k 


2 

2 


| A |  <  I.96  x  10'11 

|A| 


(2-5-18) 


After  A  has  been  determined,  we  will  be  able  to  estimate  the 
range  of  validity  for  the  approximation  (2.5-17)* 

We  will  arbitrarily  assign  n  =  0  to  the  root  to  be  obtained 
in  this  case.  Then: 


b (^1^2_^1^3^2  )  **  6 1^2  "*■  ®2^1  * 
We  approximate  y3 1  s=  k^.  Then 


(2.5-19*) 


*  63  Ckl  +  6lh(k2  '  Xo>] 


2^2 


(2. 5 -19b) 


We  assume  the  parameters  to  be 
following : 

4 

h  =  9  x  10  meters 
o^=4  mho/meter 

а,  =  lO-"3  mho/meter 
3 

6X  =  7*2  x  1010/f 

б,  =  1.8  x  105/f 

kL  =  (1-i )  3.98  x  10-5 

k?  =  (1-i)  6.28  x  10"6 
k2  =  2.09  x  10"8  f 


of  the  order  of  magnitude  of  the 


(2.5-20) 


jz 
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The  choice  of  h  and  c^  has  already  been  discussed.  The  value 
for  is  commonly  accepted  as  a  convenient  average  for  sea 
water,  although  in  certain  cases  it  might  have  to  be  revised 
downward  or  upward,  as  for  example  near  large  quantities  of 
polar  ice  in  the  first  instance,  or  in  brackish  seas  with  high 
rates  of  evaporation  in  the  second.  Any  minor  adjustments  will 
not  affect  the  validity  of  the  following  equations. 

Then : 


A 


2 

o 


h 


(2.5-21) 


From  this,  |A  |  =  - — — — —  •  For  values  of  the  parameters 

h 

1  2*5 

exactly  as  above,  this  is  given  numerically  by  — -y  .  From 
(2.5-18)  we  find; 


f  <  23^0  cps 

which  validates  (2-5-17)  well  over  the  entire  1-1000  cps  range. 
A  final  expression  for  can  be  written 

Ao  *  k2  C  1  +  t1'1)  "2?-  '  J1/2  (2.5-22a) 

1  h 

where,  of  course,  the  square  root  in  the  fourth  quadrant  is  to 
be  taken.  With  the  parameters  listed  above,  this  becomes 


X  =  k,  [  1  +  (1-i )  ]1/2  (2.5-22b) 

0  *  JI 

We  see  from  (2.5-22a)  that  the  root  approaches  ^  at  higher 
frequencies.  We  also  see  that  as  o-^  -&o  ,  XQ  ->  ^  .  Actually, 
of  course,  if  we  let  the  ionosphere  be  a  perfect  reflector  and 
the  earth  an  imperfect  one,  k  should  assume  a  value  very  close 
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to  k 2  but  with  a  small  imaginary  part  that  takes  into  account 
losses  in  the  sea.  In  our  equation  this  term  does  not  appear, 
since  it  has  been  discarded  as  negligible  in  comparison  with 
the  term  arising  from  losses  in  the  ionosphere. 

For  convenience  we  here  summarize  the  expressions  for 
Hertz  potential  in  water  due  to  the  vertical  electric  dipole: 
Branch  line  2 : 


..  2  Jl 

M  — - 7=  7 . 

6l65  b2^S  L 


iik1(z-zo)  e 


-ap  {1+1) 


(2-5-15) 


Mode : 


2iirM 


A0H0  (^P) 


in  which  |3^  =  -/k^  -  k( 


ei/5°  (z-zQ) 

<K1  +  Z2 ) '  X=X 


(2.5-23) 


2.6  The  Hertz  Potential  in  Air: 

From  equations  (2*3-1)  and  (2-3-2)  we  not  take  up  the 
branch  line  2  integral  in  a  detailed  way.  We  find 


F (ft P>2‘ ~  ^ (P ^2 '  ”^5 )  ~ 


Since 


K  K 

'  C°S  ^1+Z2^1'/32'^3)  "  K1+z2t<V|32'"V  ^ 

. 

rz2  (/3i'i32'^‘5)  Z?  1' &?' l] 

■  1  Sin^2ZlK1+Z2  ’  K1+Z2(pi,p2,-/3^)  ]  * 

(2.6-1 


K2  -  16^2  (fi^gC  -  fi-zS) _ 

Kl+Z2  ^i^55^2C"^3®  ®i^2  ^|S3C+  ^  3^  2  ^  ^ 


(2.6-2) 
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it  is  found  that  the  multiplier  of  cos  $2z  given  by 

2  ifi 


P(X) 


(2.6-5) 


where 


P(A)  =  6|)3|  (/31C+61/32S)'£  -  ^^(/31S-61/32C)4 

The  term  multiplying  sin  ^2^  is  given  by 


„  T  ^3+63^2t  ^3'65^2T  1 

’lp2  L6^2  C/5 1+6 T  )+/33 ( 6  1/32  -/&  XT )  -  6 3/32  -^3(61^2-/32T)  J 


_  2  6163^1^^3 
POO 

Expression  (2.6-1)  becomes 


(2.6-4) 


ri/3lzo  j  26i63^2^3^1sin^2z'61^2cos^2z^ 


tt2 

The  integral  is  then 

k. 


PW 


(2.6-5) 


-2M6, 6, 
t  5 


H0‘(Ap)e“1/3lZo 


k-,-ioo 

5 


j81sin£2z  ~6  2.^2cos^2z 

POO 


XdX 


(2.6-6) 

When  the  change  of  variable  X  =  a(l-i-it)  is  made  and  the 
asymptotic  form  of  the  Hankel  function  used,  this  becomes 


-i2M6163^T  e-ae(1+i)  a5/2  J^tR(t,z)  e^dt 

(2.6-7) 


where 


R(t,z)  =  £2  [t+2(l+i)]1/2(l-i-it)1/2  |^1sin^2z~6l^2cos^2Z 


-i8,z 

!  ^1  O 


-38- 


Then 


2{b1sin  b0z-6.b«cos  b~z) 

■>(0,2)  =  1  ;  ,2  _  1  2-  2  o-lblzo  (2.6-8) 


*>2  (b^+fi^S) 

Therefore,  at  sufficiently  large  distances,  this  branch  line 
integral  behaves  as 

•o/Tu  -ik,z  -aD(l+i) 

-\2J2  Mae  lo  e  /i_  _2_  -u  _  ,  -l. n__\ 

- ^ ^  ^  (b^sin  b^ ^ o  ^*^2 ^ os  j 


6 ,  6-,b~>s 
1)2 


The  mode  contribution  follows  immediately  from  (2*3-4) 


(2.6-9) 


'/3- 


2ttMA 

- 2.  H 

CT-,  o 


2/,  x.-i/^V  6lCOS  ^2°z  \^2°  I  sin  ^2°z  ^ 

(Aop)e  1  o  - (i^  +  ^2 ) '  ^ -  (2.6-10) 


X-V 


3.0  THE  SUBMERGED  HORIZONTAL  ELECTRIC  DIPOLE 


3.1  Introduction  of  Equations 

Although  the  extension  of  the  previous  results  to  the  case 
of  the  horizontal  electric  dipole  is  straightf orward,  it  is  con¬ 
siderably  more  involved  algebraically.  The  Hertz  vector  must 
be  written  with  two  components:  one  horizontal,  the  other 
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vertical,  in  order  that  the  boundary  conditions  be  satisfied.  ' ' 

We  will  take  the  dipole  orientation  to  be  along  the  x-axis. 

Then  the  two  components  of  the  Hertz  vector  are  II  and  II  .  The 

X  z 

field  components  are  derived  in  the  usual  manners 

i  =  k2  f  +  W-ff 

n  n  n  n 


k2 

n  = 

1,2,3 

(3.1-1) 

iT  . 

=  ,n  v  x  If 

n 

^|J.0  n 

match 

boundary  conditions 

we  will  need; 

E 

xn 

=  k2  n  +  |— 

n  xn  dx 

7.ffn 

n  =  1,2,3 

(3.1-2) 

E 

yn 

t> 

H 

n  =  1,2,3 

(3.1-3) 

k2  an 

Hxn 

n  zn 

iajp,o  ay- 

n  =  1,2,3 

(3.1-4) 

k2  an 

an 

Hyn 

n  f  xn 

icu(j,0  ^z 

zn\ 

•  &r-> 

n  =  1,2,3 

(3.1-5) 

^A.  Sommerfeld,  Partial  Differential  Equations,  (Academic 
Press  Inc.,  New  York,  1949)  P.  25b. 
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These  will  be  continuous  across  the  interface  if  we  require, 
from  (2.1-2)  and  (2-1-3) s 


’•ffn  *  '"K+l 

n=l, 

z=o 

(3-1-6) 

kn  "xn  “  kn  Dx(n+1) 

n=2. 

z=h 

(3.1-7) 

From  (2.1-4)  we  require 

^n  ^zn  ~  ^n+1  ^z(n+l) 

n=l, 

n=2, 

z=o 

z=h 

(3.1-8) 

and  from  ( 2 .1-5)  , 

1-2  SI1xn  _  2  5nx(n+l) 

n  3z  n+1  dz 

(3.1-9) 

We  will  first  develop  the  formulae  for  n  .  In  this  case, 
as  for  the  vertical  dipole,  there  is  primary  stimulation  in 
medium  1  only,  but  secondary  stimulation  in  all  media.  In  the 
same  way  as  for  equation  (2.1-18)  we  write: 


CO 

nx  =  M  J  Up)e-i/3i  o'  dX  (2-1-10) 

xl)pri)  0JTP1  ° 


n  =  M  fw^(x)ei^lZ  tt-  1  °  JQ(x)dX  (3.1-11) 

xl( sec)  J  x  ipl  ° 


so  that  the  total  Hertzian  vector  x-canponent  in  medium  1  is 
00 

nx  =  M  [e~i^l^z+zoU  W1(x)ei^l^z_zo^  ]  Jq(Xp)  dX 

1  0  1  (3-1-12) 


In  the  other  media: 


00 

n  r 


n  =  m 

X2 


V2(x)e_i^2z  +  W2(x)ei^2z]  e_i^lzo  JQ(XP)dX 


(3-1-13) 
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QO 

nx  =  M e_i^lzo  JQ(kp)dk  . 

3  0  1 

(3.1-14) 

Application  of  conditions  ( 3.1-7)  and  (3.1-9)  now 

yields: 

kl^l(Wl_1)  =  k2iS2(W2’  V2^  » 
^2{Vi2ei^2h-V2e-i^2h)  =  -  k^V^e-1^11  , 

k^(w1+i)  =  k^(w2+v2)  , 

k2(w2ei^2h+V2e"i^2h)  =  k^V^e_i^3h  . 

(3.1-15) 

From  these  we  find 

/3  U)  -  U2(k) 

Mi^>  -w  ♦  »,(»)  ' 

(3.1-16) 

where 

+i^2T-, 

U2  ^2  LjB2+ij5^T J  »  T  tan  ^2h  ' 

(5.1-17) 

and 

v2(k)  =  -  i6x  ^  (^)  , 

(3.1-18) 

W2(X)  -  -  iS,  , 

(3-1-19) 

»,  fj,fi2ei3;h 

V  (k)  =  2  i  — i- 4 -  . 

5  *3  (jS^M^c+i/y) 

(3-1-20) 

Evidently  U2  is  even  in  jS2,  hence  W1  and  therefore 

even  in  jS9 .  n  can  be  written  in  the  form 

n  are 

X1 
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2  * 


O.  --JL  r,p2cos  P2Z-  lu2sin  P2Z\  e~iplzo  ,  ,,.x_ 

nx,  ~  2M  rzj  ( - CTT - )  "Tisr1  J0Up)*d* 

2  k20  12  2  (3.1-21) 


from  which  we  see  that  it,  too,  is  even  in 

We  will  now  determine  the  form  of  n  .  From  (3.1-6)  we 

z 

note  that  at  the  boundaries  z  =  0  and  h. 


*\  anz(„+l)  3Iix(n+l)  SII*„ 
5^ -  dz  6x  ”  Sx 


for  n  =  1  and  2,  respectively.  However 
3nx  anx  A  snx 

ox  515  cTx  op 


(3.1-22) 


(3.1-23) 


Therefore  a  factor  cos  cp  must  occur  in  the  left  side  of  equation 
(3.1-22),  and  apparently  we  will  not  be  able  to  use  functions  of 
the  type  JQ(  ?^p)  e^^n2,  but  must  resort  to  those  of  index  1: 

cos  cp  J1(Ap)e±1^nz  . 


(3.1-24) 


No  primary  stimulation  is  present  in  this  case.  Then 
in  the  three  media  we  can  write: 


n  =  M  cos  cp 
Z1 


cp1(7^)ei^lzJ1(7vp)  A4dA  ,  (3-1-25) 


-^lZo  ,  2, 


00 

n  r* 


Z2  =  M  COS  Cp 


nw  =  m 
Z3 


[t2(70e_i^2z  +  cp2(A)ei^2z]j1(Ap)e-ij5^-  0  X2d* 

0 

CO 

cos  cp  ^(x)e~i^3zJ1(xp)e^l  0  *2d* 

r\  1 


i  0]*, 

^1 

(3.1-26) 


(3.1-27) 
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We  have  modified  the  coefficient  in  each  case  by  inclusion  of 
the  factor  ^e“^lzo/i/3^ .  This  will  simplify,  somewhat,  the 
determination  of  the  coefficients  from  the  boundary  conditions. 
We  must  satisfy  at  z=0  and  z=h: 


knnzn  “  kn+l  nz(n+l) 


(3.1-8) 


an  an2  an  an 

=  cos  CP  (d.p-(  "■•*•!?  -  gJi) 


5z - 5z" 


(3.1-28) 


for  n  =  1  and  2,  respectively.  Four  equations  are  obtained,  as 
follows: 

kJPl  a  k|(t2  +  *2)  » 


k2(t2e-i^2h  +  cP2ei^2h)  =  k|t5e_i^3h  > 


(3-1-29) 


1/3 


2A: 


2*2  -  ^2*2  +i^l  =?7T^(1  +i6l} 


-i(i2*2e-^2h+ip2®2e1P2h+ifl3*Je-1^h.-  2  j  ( jiffies) 


From  these  we  find 

20i[(  i-^)  ( /S  S-6  fl2c)  ( /S2C+i/S  S)  -6X<  i+6  ) 


*1  " 


( jS2C+i/35S)  (  S^C-^S)  (  ^4U2)  ( i31+X2) 


*2  =  -  7j5J  [(^i+i6l^2)cpl  + 
^2  =  7^  [(^i-i6i^2)cpi  + 


(3.1-30) 

(3-1-31) 

(3.1-32) 
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*5-S 


i/3,h 


fr  ^lS  +  6i|32c^cpi 


3P2 


2(i-6  )/3  S-, 
i5!- + 


(5.1-33) 


where  X2  =  (-ia-j^),  with  Z2  defined  as  in  (2.1-21).  We  'see 
that  cp^  is  even  in  and  It  is  readily  verified  that 
(ty 2e~^^2z  +  ^e1^2)  is  also  even  in 


2.2  The  Hertz  Potential  in  Region  1 

We  first  take  up  the  n  component,  hence  we  refer  to 
(2.1-12)  in  connection  with  (3.1-16).  (2.1-12)  is  changed  to 

a  contour  integral  in  the  complex  A-plane  by  precisely  the 
same  methods  that  were  used  in  the  case  of  the  vertical  elec¬ 
tric  dipole.  The  portion  of  the  integrand  which  multiplies 
the  Bessel  function  is  easily  seen  to  be  well  behaved  as 
|aI  ->  00  in  the  right  half  plane,  when  contours  are  taken 
between  poles  of  the  hyperbolic  tangent,  as  discussed  pre¬ 
viously.  Hence  the  integrands  again  vanish  faster  than  -j-iy 
on  the  appropriate  portions  of  the  infinite  semicircle  when 

I Im  A I  >0. 

We  are  again  concerned  with  branch  points  in  the  fourth 
quadrant  arising  from  the  points  were  /3^,  /3 2,  and  j6^  are  zero. 
The  integrand  is  even  in  30  that  once  more  the  problem 
reduces  to  integration  along  branch  cuts  for  j3^  and  0^,  plus 
the  determination  of  residues  at  my  poles  that  may  exist. 

The  branch  cut  for  (branch  line  2)  will  again  be  made 
vertically  downward  from  the  point  k^.  Then  the  integral 


concerned  is 


IBL2  = 


f  J  [w1(/51.^2^3)-W1(/3lf/52,-/33)]  ei/3l(Z-Zo)  ^H^ApJdA 

k,-ioo  1 

5  (3.2-1) 


Equation  (3.I-I6)  gives  the  form  of  W^. 

As  for  the  case  of  the  vertical  dipole,  the  branch  line  3 
integration  as  it  occurs  in  our  problem  can  be  expected  to  be 
negligible  at  appreciable  distances?  consequently  we  will  not 
need  a  formal  expression  for  it. 

An  expression  for  mode  contributions  due  to  the  n  component 
of  the  Hertz  vector  will  also  be  omitted.  It  will  be  found  later 
that  the  mode  represented  by  the  pole  Aq  is  not  excited  by  ..  _s 
component,  and  as  argued  previously  and  in  Appendix  A,  no  other 
modes  of  importance  are  excited. 

We  next  take  up  the  IIX  component  of  the  Hertz  potential. 

The  transformation  to  a  contour  integral  in  the  complex  A-plar.e 
proceeds  as  before,  except  that  now  we  deal  with  Hankel  functions 
of  order  one.  We  note  that 


3F  H01(z)  -  -  H^z)  , 

(3-2-2) 

3z  H02(z)  =  '  h12(z)  ' 

We  showed  in  the  analysis  for  the  vertical  electric  dipole  that 


p  ioo 


P  -i« 


H01(AP)  F (A )  AdA  =  i 


(Ap)  F(A)  AdA  (3-2-3) 


o  o 

where  F(A)  is  any  function  which  is  even  in  A.  We  take  the 
partial  derivative  with  respect  to  p,  obtaining: 
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ni00  !  9 

Hl  (Ap  )  F  (X )  X  dX  = 

o 


p  -i“»  9  9 

H^(Xp)  F(X)X*dX  , 

°  (5-2-4) 


showing  that  the  integrals  along  the  imaginary  axis  cancel  in 
this  case,  also.  Since  the  behavior  of  these  Hankel  functions 
for  large  argument  is 


(5-2-5) 


we  can  apply  the  same  reasoning  as  before  to  the  analysis  of  the 
other  portions  of  the  contour  at  radius  R,  as  R  »oo  . 

The  integrand  of  (3*1-25)  is  single  valued  around  branch 
line  1,  so  the  integration  around  branch  line  1  is  zero,  as 
before. 

The  integral  around  branch  cut  2  is  readily  obtained  by  ref¬ 
erence  to  (3.1-25).  The  expression  will  be,  after  we  have  con¬ 
verted  to  the  Hankel  function  representation: 

k^  ^  z 

iBl2  =  f  cos  0  J  [01(^1^2^3)"!?1^1^2'^3)]ei^lZHl2^P)£_ir^£A2dA* 

k,-ioo  1 

5  (5-2-6) 

The  remaining  item  of  interest  is  the  mode  solution. 

This  can  be  expressed  by: 


2tt  M  cos 
where 


h  (ft  ]_i  $2  '^3 )  ~ 


AnHl2(AnP)e^in(Z“Zo)  [fc(/3r£2,03)]  X=X  (5.2-7) 

y  n 

(i-61)(^3S-5^2C)(^2C+i^3S)  -  61(i+63)^22 

-  Tn  h  L(^2c+i^3s)(63^2c^3s)^l+u2)('3l+x2)] 
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3.3  The  Hertz  Potential  in  Region  2. 

Taking  up  first  the  nx  component  we  note  from  (3.I-I3), 
(3.I-I6),  and  (3.I-I9)  that  the  part  of  the  integrand  multiplying 
the  Bessel  function  is  well  behaved  as  |  A  |  -»  00  ,  so  that  the  con¬ 
tour  integration  can  be  performed  in  the  same  way  as  before. 

The  branch  line  2  expression  is  written: 


M 


BL2  ~  2  J 


k,-i  00 

3 


5|e-^2Z[v2(i31,i32^?)-V201^2,-^3)]  + 


(3-3-1) 


5i^2Z[w2(j31</32,J35)-W2(i31^2/-^3)]|  e_i^lZo  HQ2(Ap)d*. 


Further  evaluation  will  be  done  later. 

As  mentioned  in  the  last  section,  the  important  mode  is  not 
excited  by  the  IIX  component,  so  an  expression  for  contributions 
due  to  poles  will  be  omitted. 

In  the  case  of  the  Uz  component  we  note  from  (3.1-26), 
(3.1-3l)»  and  3*1-32)/  the  fact  that  0^(X)  is  well  behaved  as 
|^|  -*  00  ,  so  that  our  contour  integration  can  proceed  as  for 
the  region  1  discussion. 

The  form  of  the  branch  line  integral  for  II  can  be  obtained 

z 

from  (3.1-26).  Specifically: 


n 


=  ^  cos  0 


(  *2e"1/32z+02el/32z  |  (Ap)  — J^-1—  X2dX 


C  (3-3-2) 

where  C  is  the  contour  around  the  branch  cut  and  f2  and  (?2  are 
specified  by  (3. 1-31)  and(3-l-32).  We  will  leave  it  in  this 
form  for  the  present  to  avoid  writing  a  cumbersome  expression. 


Also  from  (3.1-26)  we  find  the  expression  for  the  residues 


at  the  poles.  This  is: 

_i£  nz  x 

**  *i  H12<V>  ‘-t-V5  -rn  ,  (3-5-3) 

i^l  x=xn 

where 

I  x  (X  )+n(k )  1/3 ,  sin  /3»z  -  6,n(X)/3,cos  fir>z 
yO^/Jj./^.z)  =  ± - i - ± - * - L- 

-  x- |y[  (^2c+i^3S  ^63^2C-^3S  ^  ^1+U2  ^  ^1+X2  ' 

in  which 

x(*)  =  2(i-61)(^2C+ij35S)(65i32C-i33S)(j31+X2)  , 

n(A)  =  2j31  [(i-61)05S-6^2C)02C+i/35S)-61(i+65)i322] 

3.4  Application  to  the  Water-Air-Ionosphere  Problem. 

We  shall  now  consider  the  problem  at  hand.  The  question  as 
to  whether  important  modes  of  propagation  are  excited  will  be 
taken  up  first,  since  it  has  a  bearing  on  the  fields  in  all  re¬ 
gions.  It  has  already  been  asserted  that  only  one  important 
mode  is  expected  to  exist,  due  to  the  physical  nature  of  the 
situation.  Therefore  we  need  to  test  the  integrands  of  the  ex¬ 
pressions  for  IIx  and  IIz  to  determine  whether  the  pole  X=kQ 
occurs . 

In  the  case  of  Ilx,  this  implies  that  we  must  look  for  the 
corresponding  zero  of  (0^+U2).  The  small  argument  approximation 
for  tan  jS2h  is  made,  and  the  resulting  equation  is 


When  this  is  solved  for  X,  it  is  found  that  the  value  obtained 
is  far  in  excess  of  that  which  could  be  justified  in  making 
the  approximation  to  the  tangent.  The  conclusion  is  that  the 
integrand  has  no  pole  in  the  vicinity  of  • 

Such  a  result  could  have  been  anticipated  on  the  grounds 
that  we  really  cannot  expect  a  horizontal  component  of  the 
Hertz  vector  to  excite  a  vertically  polarized  mode  of  propagation. 
II  can  give  rise  only  to  horizontal  components  of  electric  field, 
and  these  are  highly  damped  in  the  "waveguide"  region  because  of 
the  requirement  that  tangential  E-fields  be  quite  small  on  the 
interfaces . 

For  the  case  of  the  II  component,  however,  a  factor 

=  -  i  c^  (K^+Z^)  occurs  in  the  denominator.  Obviously 
the  mode  connected  with  the  pole  will  be  excited  by  this 
component. 

The  expressions  for  the  mode  contributions  in  regions  1 
and  2  will  be  written  later.  The  branch  line  integrals  will 
now  be  investigated. 

5-5  The  Hertz  Potential  in  Water 

We  take  up  the  x-component  first.  Referring  to  (5.2-1) 
it  is  found  that  the  part  of  the  integrand  in  square  brackets 
is 

_  . 

'322(^lC+ij32S  )2-'332(is2c'*'±^is)2 

It  will  again  be  necessary  to  resort  to  the  asymptotic  approxi¬ 
mation  to  the  Hankel  function 
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H 


>2(xp)  *  ^ 


o  '-  r  '  -  VirXp 
The  integral  becomes 


(Xp  -tt/4) 


(3*5-1) 


iBL2  m  2m 


ei>/4  J3  X(X)  e"i?vpdX  ,  (3*5-2) 


5 


where 


0 ,2^,X1/2e^l(z-zo) 

X(X)  =  — ^ - 5 - *—  « 

P2  (]31C+i^2S  )  -$-£  (/32C+i/31S  T 

When  the  substitution  X  =  a(l-i-it)  is  made,  (with  a  =  Re  k, 
=  Ira  kj)  the  factor  becomes 

^t  yrTHTri+rj- 

We  can  then  write  the  integral  in  the  form: 


lBL2 


2a5^2M e1^^  e-ap(1+1)  v't  G(t)e-aptdt 


(3-5-3) 


where 


G  (t ) 


Vt+2(l+i)(l-i-it)1'/2/322e1/3l^z_zo^ 
|322  (01C+i^2S)2  -  /352(^2C+iP1S)2  ' 


in  which  the  /3‘s  are  understood  to  be  expressed  as  functions  of 
t.  The  first  term  in  the  Taylor  series  expansion  of  G(t)  around 
the  point  t  =  0  is 


0(0)  = 


2  eibl<z-zo> 
(bjC'  +  ib2s)2 


where  b  =  ./k  2  +2ia2  ,  n  =  1,2;  C  =  cos  b0h,  S 
n  v  n  c. 

The  corresponding  first  term  of  the  integration  is 


sin  b2h. 


ibl2  *  ■  ^  Mel5ir/4  a  G(°)  -  ap-pr— ■  •  (5.5-4) 

-*  h 

Succeeding  terms  will  go  as  1/p  ,  1/p  ,  etc.  We  conclude  that 
the  branch  line  integral  associated  with  j3^  vanishes  at  least 
as  eap/p2  at  large  distances. 

We  will  ignore  the  IBLj  contribution,  since  it  will  con¬ 
tain  the  term  e-1*lp,  and  this  must  be  negligible  at  appreciable 
distances . 

Next  we  take  up  the  branch  line  contribution  for  Ilz  com¬ 
ponent  of  the  Hertz  potential.  We  refer  to  (3.2-6)  and  (3. 1-30) . 
Let  be  defined  by 

^1^1,^2’^3^  “  ^].(  ~  ^1  "  ^1  ’  (3.5-5) 


where  0^  is  given  by  (3.1-30)  •  Then 

0X  -  m  2  ]3161[(f-5)  -  63iS22(g-g)]  , 


(3-5-6) 


where  we  have  put  (6^-i)  t*  6^  and  (i  +  5^)  «  ^  - 
The  term  f-f  is  defined  by 


f-5  -  f(£^,£2,Aj)  “  f(^2’^2,-^3^  *  (3.5-7) 

where 

=  ( ( J 5'1'+Xp"  ’ 

with  U2  and  X2  defined  by  (3.I-I7)  and  (3.I-3O) .  The  term 
(g-g)  is  defined  by 

g-g  =  9(P1,&2>Py>  -  (3.5-8) 


where 
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It  is  found  that  to  terms  of  order 


Cf-f] 


X«k. 


65(^1+ii32T)2  (A1+61A2T)Z 


2A3(ltT2)  [jS1(61+63)  +  61A2T(63ti)] 

(5-5-9) 

7  (5.5-10) 


2^2/33G(X) 

[g-g]x*k,  -  :  2fl 


k5  53^2H01C+i^2S)4(6li32S+^1C)< 


where 


G  (M  =  (61+63)jS1^2C2  +  61^2SC(i+63)+(i63-l)^12SC+i(6163-l)|S1i32S2. 


When  approximations  based  on  the  relative  values  of  the  various 
quantities  in  the  vicinity  of  X  =  are  made,  it  is  found 
further  that 


2  )A-i  3 

-  -:j"V  2^-  +  0  033)  , 

6163^1  02  s  5 


(5.5-11) 


[g-g] 


-2  i/3, 


616^2^iC2 


+  0  035)  . 


(5.5-12) 


Terms  of  order  jSv5  are  dropped  in  our  formulation  because  of  the 
fact  that  the  Taylor  series  expansion  of  part  of  the  integrand 

is  around  the  point  $3  =  0.  Terms  of  order  /3 ^  contribute  to 

-2  7>  -3 

the  p  term;  those  of  order  contribute  to  the  p  terms, 

-2 

and  so  forth.  We  are  confining  our  attention  to  the  p  term. 
Then,  4i 


O1-0J 


X*k3  /3 2  ^ Cz 


(5.5-15) 


Our  approximations  have  reduced  this  quantity  to  quite  a  simple 
expression.  Without  such  approximations,  it  would  be  extremely 
long  and  unwieldy. 
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The  calculation  of  the  branch  line  contribution  to  II _ 


now  proceeds 


Kbl2  -  I  cos  0  J 


(of  )  e^l^z~zo^  2 

1  1  — “T75 -  X  dX  , 

C  ^1 

(3-5-14) 


where  C  denotes  a  small  part  of  the  branch  cut  in  the  vicinity 


of  k_ .  Then 
3 


2  JT  Mcos0  ei7r/4a7/2e-ap  (l+i) 


k^  JWp 


Jt  J(t)e"aptdt  , 

(3-5-15) 


in  which 


J(t)  -  ^■t+2,(,1+i)I1/yll-i-^.)3/2  e^l(z-2o)  . 

|32  <T 


Then 


J(0)  =  2U~^-  elbl^Z"Zo^  , 
b2  CT 


and  the  first  term  approximation  to  the  integral  is 


KBL2  ~  M  »  7^4 

JC,  D0  t 


ik,  Ci 


As  already  mentioned,  the  II  component  does  not  give  rise  to 

any  freely  propagating  modes,  but  the  n  component  excites 

z 

exactly  the  same  mode  as  the  vertical  electric  dipole.  The 
denominator  terms  other  than  (/3^+U2)  and  (y3 1+X2 )  merely  cancel 
the  common  denominator  of  (/3^+U2)  and  (j3^+X2).  The  determin¬ 
ation  of  the  mode  solutions  in  the  sea  then  resolves  to  finding 
the  residue  of  the  integrand  of  (3- 1-25)  at  the  pole  X=^0- 
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Referring  to  (5. 2-7)/  the  contribution  from  this  pole  is 
given  by: 


2irM  cos  0  A  H,2(A_p)  e1|3l°(z"zo^  ,r^°\  - - , _ 

0  10 

where 

_  (i-61)05S-6^2c)(/32C+ij35s)  -6  x  (i+S^2 

(6^2C-^3S)[/3102C+i^3s)  +  /32(/53  C+i^2s)‘ 


(3.5-17) 


The  factor  -  1AQ  has  been  inserted  in  the  denominator  of 

(3*5-17)  in  order  to  cancel  the  common  factor  -X  which  will 

o 

result  from  the  indicated  differentiation  and  evaluation  at 
Aq.  This  completes  our  investigation  of  the  Hertz  potential 
in  the  water. 


5.6  The  Hertz  Potential  in  Air. 

We  refer  to  equation  (3. 5-1).  Under  the  conditions  of  this 
problem,  we  find  that 


v2  (P]_«^2'^3^ -V2  (^i'^2 '  “^3) 


and  that 


2^22^3^1~^2) _ 

j322  (^C+i^S 7 -j352  (^cTi^s  )2  ' 

(3-6-1) 


w2  (^i'^2'^3^”w2  (^i'^2'  “^3)  -  — 2” 


-2/32^/3  ,(/31+^2) 


(7  • 


^2  (^ic+i^2S '  ^3  (^2C+i^ Is 

(3-6-2) 

The  branch  line  integral  is  therefore  (asymptotically) 


(3.6-3) 
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where 


X^P^/S-OjCos  ^2z+ij32  sin  £2z)  e_1^lzo 

Y  (  X  )  =  - m—  —  —  1  n - n - 15 -  . 

^2  0LC+i/32S  -  p^ip^+ip^s)* 

With  the  change  of  integration  parameter  as  before,  we  have 


lBL2  ~  2a2fl  a/fF  elT/4  e_ap(1+1)  J  JZ  S (t)e-apt  dt, 

(3.6-4) 


where 


S(t)  = 


Then 


V/a(l-i-it)1^2e~'i'^lZo  P2(p^cos  |S2z  +  i/32  sin  Pr^z)  ^t+2  (1+i ) 


S(0)  = 


P22  (PjC+ip^)2  -  TJW2 C+i^sj7 


2  ^/a  e~lblZo  (bj^cos  b2z  +  ib2sin  b2z) 


— - ~~7~ 

(b^  +  ib2S)^ 


where  the  quantities  are  defined  as  before. 

The  first  term  in  the  approximation  to  the  integral  is 
therefore ; 


!bl2  *  J7  M  el7rA  a  1/2  s(°)  6 . P~2 — " 


(3-6-5) 


This  completes  our  investigation  of  the  IIX  component. 

We  next  evaluate  the  integrals  for  Uz .  We  refer  to  (5*3-2) 
in  connection  with  (5. 1-31)  and  (5-1-32).  Since  we  have  already 
found  0(pi,p2,p^)-d!^(p^,p2,  -Pj),  the  most  involved  part  of  the 
algebra  is  done.  Making  approximations  as  before  based  on  the 
physical  value  of  the  parameters  in  the  region  of  greatest  con¬ 
tribution  to  the  integral,  we  find  that  the  integrand  of  (5*3-2) 
is  given,  to  terms  of  order  j3^,  by 


416^5 

— 2~  2  "  2 
P^P2C 


(P2  cos  P2z+ip1  sin  P2z)  X2!^2  (Xp  )e"x^lzo 

(5.6-6) 


I 


I 
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with  the  integration  to  be  taken  from  k^  -  i*  to  k^.  The 
asymptotic  form  of  the  Hankel  function  is 

hi2(ap)  ~  vSf  e’1(xP“37r/4>  .  (5.6-7) 

The  change  of  integration  parameter  is  again  X  =  a(l-i-it),  and 
the  integral  becomes  approximately 


6^M  cos  0  a 


7/2  e 


-a  p  ( 1+  i ) 

Jp 


*/t  P(t,z)e“ap 
o 


dt  , 

(5-6-8) 


where 

(i32cos  /32z+i^lsin  ^2z)(l-i-;>-t)5/2[t+2(l+i)]1/2e';L^lzo+0ft) 

P(t,z)  =  ————————  m 

^2  C 


Dropping  the  0(t)  term,  which  would  give  rise  to  terms  in  p  ^ 

b0cos  b-z+ib,  sin  b-z  ..  „ 

Fl(0,z)  =  2  (l-i)elblzo  . 

blVC 


The  first  term  asymptotic  approximation  to  the  integral  is  there¬ 
fore 


-4m  cos  06^a2e-^7r'/,^e  *blzo 


b2COs 


b2Z+ib^sin  b2Z 


b^V2" 


e-ap  (1+i) 

- 2 - 

p 


(5-6-9) 

For  the  residue  at  the  pole  we  consider  equation  (5*5-5). 
Noting  that  X2  =  at  the  pole  A=Ao,  we  obtain  for  the  residue 
at  this  pole: 


o  -i8,°z  r(A  ,  z) 

2  iirM  cos  0  A0H12(X0p)  - - — - 2 - 


Ox+x2) 


a=a. 


(5.6-10) 
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where  r(^,z)  = 

^  ^ -6^2^ )  ^  1  )^2  J  (ff^sin^2z-^^2cos  &2Z^ 

(-lAo)02c+iA5s)(«^2c"^3S^^l+u2) 

The  fields  in  air  due  to  the  submerged  horizontal  electric 
dipole  are  now  specified  by  (3*6-5)»  (5*6-9)  and  (3. 6-11). 


4.0  DERIVATION  OF  FIELD  COMPONENTS 


4.1  The  Vertical  Electric  Dipole. 

From  (2.5-23)  the  mode  solution  in  water  is 

2i7rM  i  T7  2  \  e^l  <z-2o) 

o  AoHo  ^Aop '  Tk. +Z~ ) ' " 

1  '  1  2 '  X^A 

o 

The  field  components  are  derived  as  indicated  in  (2.1-22) 
and  are  found  to  be  as  follows.  The  superscript  (v,m)  will 
be  used  to  denote  "vertical  dipole,  mode  solution",  and  (v,b) 
for  "vertical  dipole,  branch  line  solution".  In  most  cases  the 

approximation  /3^°  «  k.^  will  be  made  except  when  the  term 

2  02  2 
ki  -  /3 1  arises.  In  this  case  the  exact  difference,  X  ,  must 

of  course  be  used. 

„  (v,m)  _  2ttM  2,,  „  2^  s  elkl(z'2o) 

Eip  -  ^7  Ao  kiHi  <\>p)  w^r  x=x  ' 

o 

„  (v,m)  _  2i7TM  ,  3tt  2^  x  elkl(z~zo) 

\>*o  (V>  (*xa2)'^  • 

o 

(v,m)  _  0.„„  A  2„  2,.  x  eikl ^z-zo^ 

2  Xo  H1  (Aop^  (Kx+Z2) 


‘lz 


(4.1-1) 


H 


10 


The  branch  line  contributions  are  obtained  by  differenti¬ 
ation  of  the  expression  (2.5-15)-  The  result  is  that  the  fields 
are  given  by  appropriate  factors  multiplying  I-qi.2  *  We  find: 

Elp (V,k)  «  ak1(l-i)TBL2  , 

(4.1-2) 


e  (v,b)  rv  2  b  2\r 
Elz  *  (kl  _bi  )jbl2  * 

H  (v'b) 

H10  *  ■  a*. 


a(l-l)  ibl2 
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Therefore 


E 


H 


(v,b) 

4Mkxa^ 

lp 

M  +  - J—J 

6l63b2  s 

(v,b) 

4a/7  Ma^e"17^ 

lz 

^  7 _ 7 

616?  b2s 

(v,b) 

2MoAi  e  ^ 

0  0 

10 

« - T~2T~  ' 

b7 

ik,  (z-z  )  e'aP(l+i) 
i  1'  o'  - » -  , 


~~T 

P 


(4.1-3) 


The  fields  expressed  by  equations  (4.1-3)  must  be  added  to 
those  of  equations  (4-1-1). 

For  the  fields  in  air  we  refer  to  equations  (2.6-10)  and 
(2.6-9)  for  the  mode  solution  and  branch  line  solution,  respec¬ 
tively.  The  mode  solution  fields  become: 

o_ 


„  (v,m)  _  2ti.iT  ,  „  2/,  x_-ij3,°z_  q(Ao,Z) 

e2 p  -  “  ^  AoHi  ( V  1  0  (-1A0)(k£+z2J 


>\=*V 


(v.m)  =  _  2Mtt  a  2h  2(  )e-ipi°z( 

^  2  rT-  n  n  '  • 


o  o  '  cr 


p(^o»z) 

(-i/ao)(k1+z2)-  ' 


(4.1-4) 


„  (v,m)  _  2iM7T  ..  2  /-v  \  -ij3,°z  p(a0»z) 

H20  - 5]”  AoHl  (AoP)e  1  0  (-1A0)(K1+Z2)\=7 

where 


q(X,z)  =  6^2  sin  ^z  -  |S^  cos 
p(A,z)  =  6^  cos  p>2z  +  s^n  ^2Z  * 

The  branch  line  contribution  to  the  fields  in  air  can  be 
obtained  from  the  appropriate  differentiation  of  (2.6-9).  This 
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results  in  some  simple  modifications,  thus: 


„  (v.b)  l*,2.+i>'ViVo  e-a»<1+1> 

E9  '  ’  ’  m - y— i - 7 -  (b.cos  b,z+6..b0sin  b0z)  , 

2p  6i63b2  ^  C>  1  2  12  2 


E 


2z 


(v,b)  W7  MaVlklzo  e'3^1*1) 


6l63b2^7' 


2 -  (b^sin  b2Z-61b2COS  b2z)  , 


H 


20 


(v.b)  __ik,z_ 


(4.1-5) 


OHJ, 


TTTbTi2" 

o  l  _p  z 


1  o  - 2 -  (b^in  b2z-61b2Cos  b2z)  . 


4.2  The  Horizontal  Electric  Dipole. 

Equation  (3.5-17)  gives  the  mode  solution  in  water  for  the 
horizontal  dipole.  Equations  (2.1-22)  prescribe  the  determina¬ 
tion  of  the  fields,  which  are  found  to  be: 

rUJ 


-  MiklTco8^2[H2(  V)_H2(  v)]ei*i(2-z0)  {_1AJ{^ 


■vqr 


=  2M7TCOS0  H2(V)xo5eikl(z-zo) 


ruc) 


(-iA0)(^+x2)  • 


k=ko 

(4.1-6) 

r(Aj 


H10(h,m)  =  -^icos0  >4K(Aop)-H2(V)]eikl(2-zo)  (.-1Ao)(g1 


vqT1 


The  superscript  (h,m)  denotes  "horizontal  dipole,  mode 
solution,"  while  (h,b)  denotes  "horizontal  dipole,  branch  line 
solution." 

Contributions  from  integration  along  the  branch  cuts  are 
obtained  from  the  following  equations,  where  Ilx^  is  (3.5-4) 
and  nzlb  is  (3.5-16): 
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2 

=  cosO 


*Vb 


E  (h'b) 

ip  _  _  ap_ 

s2n  b 

-  (h,b)  _  .....  «,  /,.  2,  2  a  ^  „  b  ,  1  *  uzl 

E10  -  -  sin  0  (kL  +  -  j-  )  nxi  +  -  ■  yz8g- 


)  n  bl  +  °  “zl 
^pz  '  xl  J  3p  dz 


>2n  b  >2n  1 

(h,b)  _  y  2n  b  ^  _  3  nxl  .  3  nzl 

Elz  ’  kl  nzl  +  cos  0  ~lp8z~  +  — 


(4.1-7) 


*nzlb 


fT  (h,b)  kl2  ^  1  anzlb  j.  *  Snxlb  > 

*V  =  lajr;  C  f -50— +  sin  J 

H10<h'b)  -  -  «  ^  -  Tj 

«u(h-b)-..i..v 

It  is  readily  found 


3p 

2 

that  to  terms  of  order  p  s 
Eip(h,b)  =  cos  0  (ki2  +  2ia2)  nxib+  (i-1^!3  nzlb 

Ei0(h'b)  -  -  -in  a  kx2  nxlb  , 

Elz(h'b)  =  cos  0  kla  (1-i)  nxlb  ' 

Hlp(h'b)  '^;sin  8"I*ib  • 

H10(h'b)  =  iklcos  0  nxlb  +  (1+i)a  nzlb  ' 

Hlz(h'b)  =  sin  0  (l+i)  a  Ilxlb  . 

I  b  and  II  ,b  are  relisted  here  for  conveniences 
xl  z  l 

b  r~  iJrr/H  eikl^Z_Zo^  e'ap(l+i^ 

II  ,  «  -2  J2  a  M  elj57r/4  i - - - 7 -  , 

xl  (b^+ibgff)0  p0 

n  b  4a 2M  cos  0  ^ik,  (z-z  )  e-ap(l+i) 
nzi  *  e  1  °  — 2 — 

k  ^2  C  p 


(4.1-8) 


(4.1-9) 

(4.1-10) 
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For  the  fields  in  air  we  deal  with  (3. 6-10)  for  the  mode 
solution,  which  gives : 

E2p(h'B)-  ilM  oos  0  >'02[h02(Xop)-h22(Xop)]  e-iklzo  (^V‘?  . 

*-\> 

E22<h-">.  2iTM  coS  e  vV <V yr^r-  ■ 


x=x. 


H 

20 


(4.1-11) 

cue  ttM  cos  0  A  2[*H  2(X  p  )-H.2(X  p)|i^H 

0  o  L  o  v  oH  '  2  v  oM,J  ^  o  (j ’ 


A-A 


in  which 


(X,z) 


_[(i-61)03S-6^2c)02C+ij33S)-61(i+63)j322](picO8j32z+61P28inp22) 

( -1/A0)  (A2c+i^5S )  (6^2C~^3S  >  0*i+U2  ) 


and 


^  _)_[^i~5l)^3S~63^2C)^2C't'i/S3S)~6l^i+63)iS22]^lsinP2z~6lP2cos'32z) 

(-1Aq)  02C+i/35S )  (6^2C-i83S )  (^L+U2 ) 

For  the  branch  line  solutions  in  air  we  refer  to  (3.6-5)  and 
(3-6-9),  which  are  rewritten  here  for  convenience: 


Id 

.  e“aP  (1+i) 

(4.1-12) 

x2 

=  Ai (b^cos  b2z+xb2sin  b2z)~ - 2 -  , 

p 

b 

.  .  .e^pC^+i) 

z2 

=  A2cos  0  (b2cos  b2z+ib^sin  b2z)- - -j - 

p 

,(4.1-13) 

2J7  M  ei7r/4a  e"itolzo 

1  = 

—n— - - 

(bjC  +  ib2s) 
iiM  6La2  e“iir/4  e“iblzo 

7  = 

- v  2.  2F  2 -  ' 

t>2  C 

where 
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The  den,  vation  Qf  field  components  proceeds  similarly  to  (4.1-7)i 
giving j  to  terms  of  order  p  : 
p  (h/  -ap(l+i) 

*2p  /■  '  =  cos  0  - 2 -  E1(b1cos  b2z  +  ib?sin  bjz) 

j:  p 

(where  R^  =  (K^'  +2ia^ )  A^  4  (l+i)  a  b2  )  ' 

^  —a  p  ( 1+  i ) 

*20  |  =  -  sin  0  — — -j - k2Al(kicos  b2Z+ib2  sin  b2z) 

E271  'b) 


-ap(if-i) 

~  cos  0  - 2 - B2  ^cos  b2?.  •*  ib^  cin  b2z) 


(4.1-14) 


(where  B2  ~  +  *2  Al^  ' 


2p 


(4  „ 


=  -  sin  0  - 2 

P 


~ap  (1+i) 


-ap(l+i)  k2  b2AL 


l^U, 


-  (b^sin  b2z-ib2  cos  b2z  )  , 


*20  4  ^  =  cos  0  - - ^ 


S^^cos  b2z  +  ib1  sin  b2z) 


H, 


'where  =  a  (l+i  )A2+ib2A^ )  , 

\  -a  p  ( 1+  .i ) 

sin  0  - 2~'  —  a (l+i)  A^(i5,cus  b2z+ib2  sin  r^/. ) 

P 


i,  as  previously  noted,  b^  can  be  replaced  by  k^ 

with  negligi  , 

3  >le  error. 


5 • 0  NUMERICAL  COMPARISONS 


Although  the  field  equations  thus  far  obtained  do  not  appear 
particularly  formidable,  some  of  them  involve  rather  long  and 
tedious  calculations  in  obtaining  numerical  results.  Therefore, 
a  set  of  simplified  formulae  were  derived  from  the  field  equations 
for  the  purpose  of  easily  obtaining  numerical  results,  and  are 
presented  in  Appendix  B.  The  results  are  presented  in  the  fol¬ 
lowing  tables.  For  detailed  comparisons  between  the  different 
antennas  at  different  frequencies,  it  would  be  advisable  to  do  a 
sufficiently  complete  series  of  computations  so  that  curves  could 
be  drawn  giving  field  strengths  versus  frequency. 

The  computations  have  been  made  for  fields  in  water  at  fre¬ 
quencies  1,  10,  100,  and  1000  cps.  In  dealing  with  the  mode 
solutions  it  is  necessary  to  pick  a  distance  at  which  to  make 
comparisons,  and  this  has  been  chosen  as  1000  km.  Unfortunately, 
the  small -argument  approximations  (first  term  of  the  power  series) 
to  the  Hankel  functions  will  be  fairly  good  only  for  1  cps,  and 
the  asymptotic  approximation  only  for  1000  cps.  Here  we  have 
extended  the  small  argument  approximation  upward  to  10  cps  and 
the  asymptotic  approximation  downward  to  100  cps,  but  obviously 

the  results  cannot  be  expected  to  be  very  good  for  these  two 

2 

frequencies.  They  will  not  be  too  bad  when  only  the  HQ  function 

2  2 

is  involved,  but  will  be  somewhat  worse  for  the  and  Hg 
functions . 

We  will  list  here  numerical  values  for  A  at  the  fre- 

o 

quencies  mentioned,  since  these  will  be  used  considerably 
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in  succeeding  calculations.  In  all  cases  the  pertinent 
constants  are  as  given  by  (2-5-20). 


f /  cps 

Mod 

Xo 

Arg  Ao 

1 

3.02  X 

10-8 

-12.55° 

10 

2.4  x 

h- 

1 

O 

1— 1 

-  6.13° 

100 

2.19  x 

10-6 

-  2.12° 

1000 

2.13  x 

10-5 

-  0.78° 

that  when  p  = 

=  10°  meters. 

<  3.02  X 

10  ~2 

exp 

(-112.55°)  ' 

0.24 

exp 

(-i  6.13°) 

Ap  =  < 

0^ 

2.185 

exp 

(-i  2.12°) 

.  21.25 

exp 

(-i  0.78°)  . 

5-1  Vertical  dipole.  Mode' Solutions . 

With  the  understanding  that  in  the  following  groups  of 
four,  the  values  for  f  =  1,  10,  100,  and  1000  cps  are  from  top 
to  bottom  and  that  the  first  term  of  the  power  series  for  the 
Hahkel  function  is  used  at  f  =  1,  10  cps,  and  the  first  term 
of  the  asymptotic  expansion  is  used  at  f  *  100,  1000  cps,  the 
vertical  dipole,  mode  solution  results  are  for  2  =  0  and 
p  =  1000  km: 
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5.2  Horizontal  dipole.  Mode  solutions. 

For  the  horizontal  electric  dipole  mode  solutions  we  find 
for  the  magnitudes  at  z  =  0,  p  =  1000  km,  and  the  previously 
mentioned  approximations  to  the  Hankel  functions: 
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It  is  seen  that  for  z  *  0,  the  ratio  of  mode  fields 
from  the  horizontal  dipole  to  those  from  the  vertical  dipole 
ranges  from  the  order  of  107  to  1(P,  depending  on  frequency 
and  component.  Evidently  the  horizontal  dipole  is  more  efficient 
in  exciting  the  mode  concerned  in  this  case. 


50  Vertical  dipole,  branch  line  solutions. 

In  calculating  the  branch  line  integrals,  the  following 
will  be  useful: 
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At  p  =  1000  km,  the  magnitude  of  the  factor  — is,  for  the 
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The  field  components  at  p  =  1000  km  become,  listing  only 
magnitudes : 
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These  fields  are  negligible  compared  with  the  mode  solution 
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5.4  Horizontal  dipole.  Branch  line  solutions. 

Branch  line  contributions  for  the  horizontal  electric  dipole 
are  calculated  from  equations  (4.1-8),  and  are  determined  to  be: 
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The  branch  line  solutions  to  the  horizontal  electric  dipole 
are  evidently  much  greater  in  magnitude  than  those  for  the  ver¬ 
tical  electric  dipole,  at  low  frequencies,  and  can  therefore  be 
expected  to  have  considerable  importance  at  close  ranges.  Un¬ 
fortunately  the  approximations  made  in  evaluating  the  branch  line 
integrals  are  not  very  good  at  close  ranges.  For  validity  within 
an  order  of  magnitude,  we  might  put  bounds  as  follows  on  the 
distances : 

f  =  1  o  >  10°  meters  , 

f  =  10  p  >  3  x  10^  meters  , 

f  =  100  p  >  10^  meters, 

4 

f  =  1000  p  >  3  x  10  meters  . 


At  p  =  1000  km  the  field  components  due  to  the  branch  line 
solutions  to  the  horizontal  dipole  are,  in  magnitude: 
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6.0  ANTENNAS 


6.1  General  Discussion. 

The  requirements  for  a  submerged  ELF  electric  antenna  can 
be  summed  up  qualitatively  in  a  rather  brief  manner.  The  funda¬ 
mental  necessity  is  that  of  obtaining  the  greatest  possible 
electric  moment.  This  moment  is  proportional  to  the  product  of 
current  and  antenna  length;  hence  the  goal  is  to  use  the  greatest 
current  and  greatest  length  possible.  Because,  however,  the 
surrounding  conducting  medium  imposes  a  penalty  in  the  form  of 
high  losses  in  any  region  where  field  intensity  is  high,  the 
maximization  of  current  must  be  achieved  by  some  method  which 
maintains  field  intensity  as  lew  as  possible  everywhere  in  the 
water. 

This  problem  can  be  viewed  from  another  standpoint,  namely 
that  of  impedance  matching.  The  intrinsic  wave  impedance  of  sea 
water  is  very  low,  which  is  essentially  the  same  as  saying  that 
the  ratio  of  electric  to  magnetic  field  is  low.  For  good  im¬ 
pedance  matching  from  antenna  to  water,  it  follows  that  the 
electric  fields  should  be  kept  low  while  at  the  same  time 
generating  currents  of  the  maximum  possible  magnitude. 

If  the  additional  requirement  is  made  that  the  antenna  be 
readily  transportable,  other  restrictions  are  immediately 
evident.  Extremely  large  spherical,  spheroidal,  or  other  shapes 
with  cross  sections  which  preclude  good  hydrodynamic  properties 
are  excluded. 
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In  view  of  the  above  considerations  it  would  seem  that  the 
only  really  feasible  submarine  antenna  for  ELF  is  one  involving 
a  trailing  wire  of  some  sort.  Under  various  circumstances,  ob¬ 
jections  can  be  raised  to  such  an  antenna;  nevertheless,  there 
seems  to  be  no  other  type  for  which  much  more  serious  disad¬ 
vantages  do  not  appear.  In  the  following,  we  will  confine  our 
attentions  to  this  kind  of  antenna  and  more  specifically  to  the 
"coaxial"  antenna  about  to  be  described. 

The  problem  has  already  been  considered  in  considerable  de- 

4 

taxi  by  Moore,  who  introduced  the  coaxial  antenna.  This  is  an 
insulated  length  of  wire,  the  input  end  of  which  is  fastened  to 
the  submarine.  The  trailing  end  may  be  either  insulated  or  ex¬ 
posed  to  the  sea  water.  In  the  first  case  the  system  resembles 
a  coaxial  transmission  line  with  the  load  end  open  circuited. 

In  the  second  case  it  resembles  a  short-circuited  coaxial  trans¬ 
mission  line.  To  obtain  an  effective  short  circuit  in  the 
latter  case,  the  amount  of  exposed  wire  necessary  at  the  trail¬ 
ing  end,  as  shown  by  Moore,  need  only  be  around  one -eighth  wave¬ 
length  (in  sea  water). 

For  any  reasonable  overall  length,  considering  the  fre¬ 
quency  range  with  which  we  are  dealing,  the  short  circuited 
antenna  will  carry  essentially  uniform  current  throughout  its 
length.  This  can  be  considered  due  to  the  high  velocity  of 
propagation  of  this  "transmission  line".  Such  an  antenna  may, 
however,  be  of  appreciable  length  compared  with  a  wavelength  in 
sea  water.  At  1000  cps,  for  example,  a  half  wavelength  in  sea 
water  is  only  about  17  meters.  At  10  cycles  per  second  it  is 
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ten  times  as  great,  or  170  meters.  Although  this  latter  figure 
might  appear  to  be  a  rather  cumbersome  length  for  a  trailing 
wire,  it  is  perhaps  not  beyond  the  realm  of  possibility. 

The  distribution  of  current  in  the  open  circuited  coaxial 
antenna  can  be  approximated  as  varying  linearly  from  its  maximum 
at  the  submarine  or  input  end  to  zero  at  the  trailing  end.  The 
effective  current  in  this  case  is  therefore  half  of  the  input 
current.  Since  the  reactance  of  the  open  circuited  antenna  is 
quite  high  and  the  difficulties  in  matching  to  it  consequently 
very  great,  it  does  not  seem  advantageous  to  discuss  it  further. 

At  constant  antenna  current  the  fields  generated  by  the 
coaxial  antenna  are  shown  by  Moore  to  become  significantly 
greater  as  the  diameter  of  the  antenna  is  increased  beyond  about 
r,  where  A^  is  wavelength  in  sea  water.  The  effect  can  be 
described  in  terms  of  an  "equivalent  current"  which  increases 
with  diameter  of  the  antenna.  For  lesser  diameters  there  is 
no  appreciable  variation  with  diameter.  In  the  ELF  range,  it 
will  be  assumed  that  cross-sectional  dimensions  equal  to  or 
greater  than  A^/tt  are  excluded  for  reasons  of  physical  bulk. 

The  equivalent  current  used  in  computing  fields  for  the  antenna 
is  then  just  the  same  as  the  actual  current  in  the  conductor. 

6.2  Power  Considerations. 

Although  the  gain  of  the  short  circuited  coaxial  antenna 
can  be  modified  by  altering  its  length  in  various  ways,  the 
effect  is  one  which  involves,  at  most,  a  few  decibels  of  power 
in  one  direction  or  another  from  the  antenna,  and  does  not 
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merit  discussion  for  our  purposes.  It  will  be  sufficient  if  we 

assume  for  any  of  the  frequencies  in  the  given  range  that  the 

antenna  is  a  half  wavelength  or  less,  so  that  the  gain  in  any 

direction  differs  from  that  of  an  isotropic  radiator  by  no  more 

than  a  factor  of  two.  We  will  also  ignore  the  directivity  gains 

which  can  be  achieved  by  the  use  of  multiple  radiating  elements, 

or  by  suppression  of  alternate  phases  on  a  line  which  is  several 

half  wavelengths.  It  is  highly  doubtful  whether  such  schemes 

can  effect  changes  of  greater  than  an  order  of  magnitude  or  so 

in  power  transmitted  between  two  points,  without  going  beyond 

the  limits  of  physical  feasibility. 

In  the  following,  then,  the  equivalent  moment  of  the  short 

circuited  coaxial  antenna  will  be  taken  as  I  l,  where  I  is  the 

o  o 

r.m.s.  input  current  and  l  the  length  of  the  insulated  section. 
(Radiation  from  the  uninsulated  section  will  be  disregarded.) 

The  impedance  of  such  an  antenna  is  shown  by  Moore  to  be 
defined  by: 


R  -  Sgi  ,  (6.2-1) 

x  =  ISr1  [°-116  "  log  +  log  b/a  ] 


=  [0.116  -  log  6 a/7  ]  .  (6.2-2) 

in  which  6  =  vc mo/2  ,  b  is  the  outer  radius  of  the  insulating 
material  of  the  antenna,  and  a  the  radius  of  the  conductor  or 
antenna  core.  The  tacit  assumption  has  been  made  that  the 
antenna  is  sufficiently  removed  from  the  boundary  between  sea 
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and  air  so  that  the  impedance  is  practically  the  same  as  that  of 
the  antenna  in  homogeneous  sea  water. 

We  will  now  consider  the  transmission  of  power  between  two 
short  circuited  coaxial  antennas.  The  power  transmitted  is  given 
by 


P 


t 


(6.2-3) 


where  It  is  the  r.m.s.  current  and  Rfc  the  radiation  resistance. 
The  received  power  in  a  receiver  matched  to  the  receiving  antenna 
is  o  o 


Pr  = 


V 

4r 


(e  iry 

4r_ 


(6.2-4) 


r  r 

where  Vr  is  the  induced  voltage  at  the  terminals  of  the  receiver, 
and  is  given  by  E^r,  with  E  the  electric  field  parallel  to  the 
axis  of  the  antenna  and  l ,  the  receiving  antenna  length.  The 
receiving  antenna  radiation  resistance  is  Rr» 

The  ratio  of  received  to  transmitted  power  is  then  given 


by 


E2  -tr2 

“rVt2 


(6.2-5) 


Since  the  field  expressions  as  previously  derived  are 
expressed  in  terms  of  M,  where 


It  <t  ,  (6.2-6) 

M  *  ~Wo~ 

it  will  be  convenient  to  write 

Pr  (E 2M2)  lr2lt2 

Pt  647r2a2RrRt 
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Using  the  expression  for  radiation  resistance  of  the  short 
circuited  coaxial  antenna,  this  is  equivalent  to 


(E2/k2)  i2i2 

W?f2 


(6.2-8) 


Letting  a  =  4, 


1.01  x  10 
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In  considering  the  fields  in  water  due  to  the  mode  propagation 
for  either  the  vertical  or  horizontal  dipole,  it  may  be  noted 
that  the  radial  component  is  the  largest.  Because  the  effect 
of  depth  can  be  computed  separately  from  that  of  radial  dis¬ 
tance,  the  receiving  antenna  can  be  assumed  to  be  located  at 
z  *  0,  regardless  of  the  fact  tha  t  the  formula  for  radiation 
resistance  is  not  then  applicable.  The  attenuation  resulting 
from  location  at  z  =  0  can  then  simply  be  added  to  that  due  to 
distance. 

With  alignment  of  the  receiving  antenna  in  the  radial 
direction,  we  find  at  1  cps  and  1000  km. 


Pr  \  (v,m) 

P t  I  1  cps 


5.1  x  10"^°  lrlt  e~2azo 


(6.2-10) 


where  a  =  imaginary  part  of  k^. 

The  height  of  the  vertical  transmitting  antenna  may  well  be  an 
important  matter,  however.  This  is  taken  care  of  by  integrating 
over  the  range  of  zQ. 
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At  1  cps  and  1000  km  we  have 
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Lo  i  i 

If  the  antenna  is  sufficiently  long  (say,  a  half  wavelength) 
so  that 

e_az2  »  e~azl  ,  we  can  write 


1.75  *  10"24  M  |- — yJL 


'“Ip1 

The  power  expression  then  becomes 

/  — |  <V'm)  =  <-r  e-2«2  , 

v  Pt  /  1  cps  k.z 

=  9-78  X  10'56  lr  e~2az2  . 

Letting  l r  be  a  half  wavelength  (540  meters),  we  obtain 
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Pfc/  1  cps 


(6.2-14) 


representing  a  system  attenuation  of  about  323  db,  not  including 
losses  due  to  submergence;  that  is,  losses  which  are  a  function 
of  z  and  %2  alone.  The  latter  can  be  calculated  readily,  and 
will  be  ignored  for  the  present. 

Let  us  compare  these  figures  for  1  cps  with  those  which 
would  obtain  at  1000  cps.  Again  with  the  antenna  assumed  long, 


Pi  (v,m) 

Pfc /  1000  cps 


2.68  x  10-54  L  e~2az2  . 


(6.2-15) 
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At  1  cps  and  1000  km  we  have 


E 


Id 


1.75  x  10 


M 


-24  ..,e'az2  -  e~azl 

I  “ 


(6.2-11) 


If  the  antenna  is  sufficiently  long  (say,  a  half  wavelength) 
so  that 

e"az2  »  e"azl  ,  we  can  write 


„~az2 

|Ei0l  -  1-75  *  1°  M  I — iq- 


The  power  expression  then  becomes 


pr)  (v'm)  3.1  x  10'40  4r  -2<xz, 

—  I  =  - - ~ - s.  e  t 

1- 


Pt  /  1  cps 


=  9.78  X  10"56  lr  e~2az2 


(6.2-12) 


(6.2-13) 


Letting  tr  be  a  half  wavelength  (540  meters),  we  obtain 


(!r|  (v'm)  .  5>2?  x  10-53  e-2az2  #  (6.2-14) 

\Ptl  1  cps 

representing  a  system  attenuation  of  about  323  r  (t  including 
losses  due  to  submergence;  that  is,  losses  which  are  a  function 
of  z  and  z2  alone.  The  latter  can  be  calculated  readily,  and 
will  be  ignored  for  the  present. 

Let  us  compare  these  figures  for  1  cps  with  those  which 
would  obtain  at  1000  cps.  Again  with  the  antenna  assumed  long. 


Pr\  (v,m) 

P fc/  1000  cps 


2.68  x  10-54  lr  e"2az2  . 


(6.2-15) 
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With  a  receiving  antenna  17  meters  (a  half  wavelength)  long, 
the  ratio  is 

4.55  x  10"-5-5  e“2az2  , 

giving  a  system  attenuation  of  about  323  db  excluding  water 
losses . 

In  the  case  of  the  mode  solutions  for  the  horizontal  dipole 
we  find,  under  the  same  assumptions: 


Pr|  (h,m) 
Pfc  1  cps 


1.22  x  10-26  il  e“2azo  , 


(6.2-16) 


which  for  half  wavelength  antennas  and  0*0  becomes 
3.56  x  10"21  e“2azo  . 

System  attenuation  excluding  water  losses  is  thus  about  204.5  db. 
At  1000  cps. 


(h,m) 

Pt  I  1000  cps 


7.56  x  10“28  lrlt  e~2azo  cos20, 


which  for  a  half  wavelength  and  0  =  0  is 
2.12  x  10“25  e'2azo  , 


(6.2-17) 


or  about  246.7  db  attenuation  excluding  water  losses. 

It  will  be  of  interest  to  compare  these  results  with  the 
branch  line  solutions.  In  all  cases  the  branch  line  solutions 
at  1000  km  have  their  greatest  magnitudes  at  1  cps,  and  become 
negligibly  small  at  1000  cps  in  comparison  with  the  mode 


solutions . 
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The  electric  field  magnitudes  are  again  greatest  in  the 
radial  direction,  although  for  the  horizontal  dipole  the  0  com¬ 
ponent  is  nearly  equal  to  the  radial.  In  any  case,  it  will 
suffice  to  consider  only  the  radial  component.  Taking  first 
the  vertical  transmitting  antenna,  but  again  with  receiving 
antenna  at  z  =  0  and  aligned  along  a  radius,  we  have 


Pr  \  (v,b) 
/  1  cps 


=  2.59  x  10-57  ^  .  e"2azo  . 

r  t 


(6.2-18) 


This  is  negligible  compared  with  the  mode  solution,  and  obviously 
the  higher  frequencies  are  of  even  less  interest. 

For  the  horizontal  transmitting  antenna  at  1  cps,  under  the 
previous  assumptions 


(h,b)  -^r.  Ortr*  9 

=  2.22  x  lO"^  i  i  e~*azo  cosz0  .  (6.2-19) 

1  cps  1 

With  lr  =  =  540  meters  and  0=0  this  becomes 

6.48  x  10"26  e_2azo 

or  about  252  db  attenuation  not  including  water  losses.  The 
horizontal  mode  solution  gives  the  least  attenuation  (205  db) 
of  any  considered,  but  it  should  be  kept  in  mind  that  more 
realistic  antenna  lengths  would  cause  this  to  be  increased 
somewhat.  With  antennas  54  meters  long,  for  example,  the 
system  attenuation  in  this  case  would  be  225  db. 

For  these  figures  to  have  practical  significance,  they 
must  be  related  to  expected  noise  power.  Unfortunately,  the 
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exact  magnitude  of  atmospheric  noise  is  not  well  known  in  the 

frequency  range  under  discussion.  In  measurements  that  have 

been  made  '  ■*'  noise  fields  have  been  given  ranging  from 

about  1  to  1800  microvolts  per  meter  per  cps  for  frequencies 

from  1  to  800  cps.  Aarons  presents  recordings  which  show 

noise  of  the  order  0.1  to  250  microvolts  per  meter  per  cps,  with 

considerable  fluctuations  about  these  values.  Willis  measured 

-8 

magnetic  fluctuations  of  6  x  10”  gauss  per  cps  at  5  cps,  ranging 
down  to  4  x  10-11  gauss  per  cps  at  800  cps.  This  gives  from 
1800  to  1.2  microvolts  per  meter  vertical  electric  field. 

_7 

Goldberg  mentions  magnetic  components  of  the  order  3  x  10 
gauss,  which  would  be  9000  microvolts  per  meter  electric  field, 
but  it  appears --without  being  entirely  clear--that  this  is  for 
his  entire  pass  band  of  1.0  to  150  cps.  If  uniform  across  the 
band  this  would  represent  noise  strength  of  60  microvolts  per 
meter  per  cps. 

Since  the  most  detailed --although  somewhat  difficult  to 
interpret — results  are  those  presented  by  Aarons,  we  will  take 
the  figure  2  microvolts/meter/cps  from  one  of  his  charts  as 
being  roughly  representative  at  frequencies  near  the  low  end 
of  the  1  -  1000  cps  range.  It  should  be  kept  in  mind  that 
this  may  be  extremely  optimistic.  We  will  consider  a  trans¬ 
mitter  of  100  kilowatts  output  power  and  a  frequency  of  1  cps, 

°J.  Aarons,  Journal  of  Geophysical  Research,  6l,  647  (1956) 
^H.  F.  Willis,  Nature,  161,  887  (19^8) 

20P.  A.  Goldberg,  Nature,  177,  1219  (1956) 
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with  radial  distance  1000  km  between  transmitter  and  receiver. 
Assuming  0=0  and  half  wave  (540  meter)  antennas,  the  received 
power  at  z  =  0  will  be,  essentially, 

Pr(h,m)  =  5>56  x  10-l6  e-2azQ  _  (6.2-20) 

When  the  vertical  electric  field  in  air  is  2  microvolts 
per  meter,  the  magnitude  of  the  tangential  magnetic  field  at 
the  water  surface  will  be 


|H 


tan 1 


'vert 1 
^o/eo 


377 


(6.2-21) 


The 


tangential  electric  field  will  be  given  by 

'  '  '  !HtaJ  ■ 


lEta„l  -  lZo' 


where  Z_  = 


(1+i)  ^ 
v  2a 


Jtan 1 


Therefore, 
1.4  x  10~3  yT 


377 


E. 


n-6 


=  3.76  x  10"°  JI  |Ev|  . 


(6.2-22) 


The  figure  2  nv/meter  at  1  cps  then  gives 

2 

|Etan|  =  7-5  x  10”  volts/meter/cps . 

Received  noise  power  with  the  antenna  aligned  for  maximum 


'  Etan '  is: 


P  = 
r 


lEta„l V  2'5  »  W5  Er\ 


4r 


(6.2-23) 


which  for  the  above  figures  becomes 

■15 


Pr  =  7*63  x  10"  3  watts. 


Thus  the  signal  to  noise  ratio  will  be 

S  |  (h,m)  _  3.56  x  10"16 

!T  1  1  cps  _  ,  -_15 

1  1000  km  7 -63  X  10 


-2az^ 
e  o 


4.66  x  10~2  e-2azo 


(6.2-24) 
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or  about  -  13*3  db.  The  number  of  decibels  of  attenuation  due 
to  the  depth  of  the  transmitting  antenna  will  of  course  subtract 
directly  from  this,  that  is,  will  be  deleterious  to  the  S/fcr  ratio. 
As  mentioned  previously,  the  figure  taken  for  noise  magnitude  is 
probably  optimistic  to  begin  with. 

It  may  be  noted  here  that  receiver  thermal  noise  is  not  im¬ 
portant  compared  with  atmospheric  noise  in  this  case.  For  example 
if  one  considers  the  receiver  as  matched  to  the  antenna,  then  for 
1  cps  bandwidth  the  noise  power  in  the  receiver  input  resistance 
is  2kT  (k  =  Boltzmann's  constant,  T  =  temperature,  degrees 
Kelvin).  For  20°  C  this  is  8.03  x  10 ~2^  watts,  which  is  about 
250  decibels  below  the  transmitted  power  of  10^  watts.  Thus 
the  basic  limitation  is  evidently  that  associated  with  natural 
noise. 

6.3  Antenna  Depth  Considerations. 

For  the  mode  solutions,  the  fields  increase  monotonically 
with  frequency,  but  with  a  dependence  on  frequency  that  is 
somewhat  involved.  Moreover,  the  dependence  is  different  for 
different  components  of  the  field,  and  varies  with  distance. 

In  Figures  6.1  and  6.2,  the  field  components  are  drawn  as 
a  function  of  frequency  on  logarithmic  graphs.  As  a  very  rough 
approximation  at  1000  km  distance,  the  dependence  might  be  taken 
to  be  (aside  from  the  exponential  term)  for  a  radial  component: 

Low  Frequency:  ~  f 7  f0’2 

(v,m)  ^  f 2 .  E  (n,m)  ^  f3/2 
lp  lp 


High  Frequency: 


(6.3-1) 
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Let  p  be  the  power  of  frequency  involved  as  a  very  rough  approx¬ 
imation.  The  optimum  frequency  for  a  given  field  component  will 
be  given  by  differentiating  the  field  component  with  respect  to  f 
and  setting  the  result  equal  to  zero.  We  thus  obtain  (assuming 
receiving  antenna  depth  is  zero) 

dk 

pfP-i  e-azo  _  iZofP  =  0  (6.3-2) 


From  this  it  is  found  that 

£  _  1.22  p  x  10 

£°pt - T2 - 


5 


(6.3-3) 


o 

For  the  (h,m)  solution  and  a  receiving  antenna  at  z  <  0, 

2  2 
zq  in  this  equation  can  be  replaced  simply  by  (z-zq)  ,  so  that 

the  optimum  frequency,  to  a  crude  approximation,  varies  in¬ 
versely  as  the  square  of  the  combined  depth  of  receiving  and 
transmitting  antennas . 

The  same  relationship  was  found  by  Moore  to  hold  in  a  more 
exact  fashion  when  the  atmosphere  is  considered  homogeneous. 

The  z-dependence  for  the  (v,m)  situation  is  more  complicated, 
but  essentially  of  exponential  behavior,  so  that  again  an  approxi¬ 
mate  idea  of  optimum  frequency  can  be  obtained  in  the  same  manner. 

As  an  example,  if  the  combined  depth  of  receiving  and  trans¬ 
mitting  antennas  is  100  meters,  one  finds  the  optimum  frequency 

is  about  12.2  p  for  a  distance  of  1000  km.  Considering  the  E. 

XP 

component,  this  would  mean  f  t  would  be  roughly  20  cps  for  the 
vertical  transmitting  antenna.  For  the  horizontal  case,  f  t 
would  be  somewhere  in  the  vicinity  of  2.5  to  5«0  cps;  and  not 
very  sharply  defined  on  account  of  the  near-zero  frequency 
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dependence.  Evidently  quite  low  frequencies  must  be  used  when 
the  antenna  depths  are  appreciable. 

The  branch  line  solutions,  on  the  other  hand,  decrease 
monotonically  with  frequency.  For  them  there  is  no  optimum, 
and  the  "best"  frequency  is  the  lowest  that  can  be  used. 

6.4  Summary  on  Antennas. 

In  the  frequency  range  1  -  1000  cps ,  the  system  attenuation 
for  a  distance  of  1000  kilometers  between  submerged  horizontal 
short-circuited  coaxial  antennas  is  large  but  does  not  preclude 
the  possibility  of  practical  communication.  The  noise  situation 
is  not  well  known,  but  at  best  it  would  appear  that  signal  will 
be  about  1J  db  or  more  below  noise.  Further  degradation  of  the 
signal-noise  ratio  will  occur  due  to  vertical  travel  of  the 
waves  between  antennas  and  water  surface.  The  quoted  1J  db  figure 
is  that  obtained  for  the  mode  solution  for  the  horizontal 
transmitting  antenna  at  a  frequency  of  1  cps. 

At  distances  less  than  1000  km  the  so-called  "branch  line" 
solution  becomes  comparable  in  magnitude  to  the  mode  solution, 
but  unfortunately  our  analysis  does  not  reveal  what  sort  of 
behavior  can  be  expected  close  to  the  antenna  for  the  "branch 
line"  solutions.  One  method  of  tackling  this  question  would  be 
by  means  of  numerical  integration,  possibly  programmed  for  an 
automatic  computer.  A  possible  approach  would  be  to  draw  a 
branch  cut  from  k,  toward  k.  so  that  Bessel  functions  of 

j  1 

argument  (e'^^x )  could  be  used.  Close  to  k^,  the  contribution 
to  the  integral  will  be  negligible  for  all  except  such  low 


-88- 


values  of  p^/f  as  to  be  essentially  without  interest.  This  is 
because  the  Hankel  function  in  the  vicinity  of  k^  becomes  pro¬ 
portional  to  e"lklp,  provided  (roughly)  pjf  >  200.  Whether 
this  approach  is  better  than  that  of  putting  the  branch  cut 
vertically  downward  would  have  to  be  determined  by  a  more 
detailed  study. 

The  branch  line  solution  is  much  smaller  than  the  mode 
solution  at  1000  km  distance  and  1  cps  and  becomes  less 
important  for  higher  frequencies  at  this  same  distance. 

The  very  slight  attenuation  with  distance  of  the  mode 
solution  would  appear  to  make  it  ideal  for  communications  pur¬ 
poses.  However,  it  requires  a  very  efficient  submerged  antenna 
to  launch  a  wave  up  through  the  sea -air  interface  before  it 
even  goes  between  the  sea  and  the  ionosphere.  It  has  been 
shown  for  a  practical  length  antenna  of  54  meters  that  possible 
communications  would  be  at  a  range  much  less  than  1000  km. 

One  does  not  exclude  a  slight  possibility  that  a  new  and 
perhaps  radical  antenna  design  might  enable  much  more  efficient 
coupling  of  electromagnetic  energy  to  the  water.  This  possi¬ 
bility  has  been  dubbed  "slight",  however,  because  the  physical 
nature  of  the  problem  seems  to  mitigate  against  it  in  every  way. 
The  requirement  of  physical  transportability  is  in  direct  con¬ 
flict  with  the  requirement  that  low  current  densities  be  main¬ 
tained  in  the  vicinity  of  the  antenna.  A  very  analogous  situa¬ 
tion  is  one  which  arises  in  acoustics,  namely  that  of  matching  the 
high  impedance  of  an  amplifier-loudspeaker  combination  to  the 
low  impedance  of  air.  The  problem  is  particularly  acute  at  long 
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wavelengths  and  is  met  by  using  one  or  more  loudspeakers 
which  are  as  large  in  area  as  is  physically  and  economically 
feasible.  Even  then,  it  is  customary  to  enhance  the  impedance 
match  by  using  baffles  or  horns  of  various  types --thereby 
increasing  even  more  the  effective  area  of  energy  coupling. 

It  is  the  opinion  of  the  writer  that  the  requirement  of 
large  size  is  just  as  basic  to  the  submerged  antenna  problem 
as  it  is  to  this. 

The  same  remarks  are  applicable  to  the  receiving  antenna. 
Because  of  the  impedances  involved,  the  efficiency  with  which 
it  intercepts  electrical  energy  in  the  water  must  depend  on 
its  physical  size. 
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When  the  distance  between  perfectly  conducting  planes  of  a 
parallel  plane  waveguide  is  less  than  a  half  wavelength,  the 
propagation  constants  for  TE  and  TM  modes  are  pure  imaginary. 

In  the  development  of  this  paper,  their  location  would  be  repre¬ 
sented  by  poles  on  the  negative  imaginary  axis  of  the  complex 
X -plane. 

To  obtain  an  estimate  as  to  the  amount  by  which  such  poles 
may  move  when  finite  conductivity  is  introduced  in  the  boundaries, 
we  will  consider  a  rather  idealized  case  where  a  wave  of  infinite 
transverse  extent  is  assumed  to  be  propagating  between  the  in¬ 
finite  parallel  planes.  The  conductivity  of  both  planes  will  be 
taken  as  equal  to  that  of  the  ionosphere;  thus  it  is  to  be  ex¬ 
pected  that  our  results  will  represent  a  worse  situation  than 
the  one  we  have  been  investigating. 

The  direction  of  propagation  will  be  taken  along  the  z- 
axis.  The  planes  are  the  surfaces  y  =  0  and  y  =  h.  For  TM 
waves  we  can  assume  a  magnetic  field  of  the  type 

Hx  =  Ho  cos  ky  y  eicut-Vz  (A-la) 

and  derive  the  remaining  field  components  by  use  of  Maxwell's 
equations.  It  is  found  that: 


*.  ■  «o  ■*»  <\  (*-“>) 

Ey  -  -  i  Ho  cos  <ky  y)eia,t-V2  (A-lc) 
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All  other  field  components  are  zero.  From  the  wave  equation 
one  obtains  the  relation 


ky2  =  y2  -  0)2/c2 


(A -2) 


If  the  walls  were  perfectly  conducting,  it  would  be  necessary 
that  ky  have  the  eigenvalues  nir/h,  n  =  1,  2,  3»  ----,  in  order 
for  the  boundary  conditions  to  be  satisfied.  In  this  case  we 
will  assume  that  the  conductivity  is  sufficiently  good  so  that 
the  eigenvalues  are  perturbed  only  by  a  small  amount  p.  Then 

ky  *  +  P  (A -3) 

and  to  terms  of  order  p, 

^2  ^  (n*)2  +  (*-4) 

For  the  fields  in  the  conducting  medium  we  can  take 

H  =  Ce-by"YZ+icut 


E  +  b/aCe-by"YZ+ia* 
z 

Ey  =  -V/aCe-by-YZ+icut 

where  satisfaction  of  the  wave  equation  requires  that 


(A -5) 


2  2 

b  +  y  -ia*i  a  =  0 


(A -6) 


Here  the  displacement  current  has  been  assumed  negligible  in 
comparison  with  the  conduction  current.  This  is  valid  for 
o/cue  »  1,  which  is  true  in  our  case. 

At  y  =  b  the  boundary  condition  requires 

sin  (k  h)ei^t-YZ  =  b  -bh-YZ+icot  (a _7a) 

H0  cos  (Xy  h)eit5t-vz  -  C  e-bh-YZ+l<»t  (A-7b) 
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Dividing  the  top  equation  by  the  bottom  one: 


ky  sin(ky  h)  iooeb 

cos  (ky  h)  o 


(A -8) 


The  mode  least  attenuated  will  be  that  for  which  n  =  1; 


hence 


k  =  J—  +  d 

V  h  * 


(£+  P)  tan  [(J+  p)  h]  -  3SS6.  . 

Approx ima  te ly , 

(ir/h  +  p)  [tan  n  +  ph  sec^  tt]  = 

from  which,  to  terms  of  order  p, 
iooeb 


P  mt 


We  can  write 


7 TCT 


b  =  J  iowc  -  (£)2  -  + 


CO 

7 


(A -9) 

(A -10) 

(A-ll) 

(A -12) 


For  frequencies  well  into  the  "stop  band",  that  is,  well  below 
the  so-called  cutoff  frequency,  the  last  two  terms  can  be 
neglected,  so 


b  V  icuua  -  (^)2 


Therefore, 

Y  M 


Jo2  -4  + 1 

c 


Ea"  v'iava  ' 


(a -15) 


(A-14) 


If  the  conductivity  is  allowed  to  become  infinite,  obviously 
one  recovers  the  expression 


=  Jo- 


/TT\2  00 

w  '7 


(A -15) 


which  is  the  usual  equation  in  such  a  case.  We  shall  now  in¬ 
quire  into  the  magnitude  of  what  me  might  call  the  "correction 


term",  C 


<s>2 


“2 1 

77 


(a -16) 

we  can  be 


roya 

If  this  is  small  in  comparison  with 
satisfied  that  a  pole  in  the  complex  plane  representing  the 
location  of  y  will  not  be  greatly  shifted  from  its  value  for 
perfect  conductivity. 

It  will  suffice  to  check  the  frequencies  f  =  1  and  1000  cps 


For  f  =  1, 


C 


o 


2cueir 

^7 


(A -17) 


and 

V  H-  (1  —  HE-)  w  H-  ^  x  10  11) 

Evidently  y  is  not  changed  very  much  in  this  case. 

CQ  *  -  3.52  x  10"7  ( 1 — i )  , 

and  one  obtains 


(A -18) 

For  f  -  1000 
(A -19) 


y  «  2.85  x  10 "5  -  (l-i)  3.52  x  10 "7  (A -20 ) 

so  that  again  the  perturbation  is  seen  to  be  slight. 

For  TE  waves  we  can  carry  out  a  similar  development. 

The  fields  in  air  can  be  written 

Ex  =  E=  sin  (ky  y)e“yz+l“t 

By-  hs  *0  ■*»  <*y  y)e-yz+iait  (A -21) 

Hz  -  %  Ec  cos  <ky  y)e-y»-‘ 


-V- 


and  those  in  the  conductor  are: 
E  =  Ae-bY-YZ+i“>t 


H  -  a  Y  -by-vz+iuit 
Hy  ”  A  HEfT  e 

H  =  ^  JL.  e-by-Yz+i^t 

Z  K41 


At  the  boundary  y  =  h. 


EQ  sin  (ky  h)  =  Ae 


-bh 


_  /.  ,  \  A  b  -bh 

E  cox  (k  h )  =  -Ai -  e 

IC41  o  v  y  '  icqi 

Dividing  the  top  equation  by  the  bottom  one: 

k 

tan  (ky  h)  =  - 

Again  letting  ky  =  ^  +  p  for  the  n  =  1  mode, 
tan  h  (J  +  p)  «  tan  ph  -  -  5-  (J-  +  p) 


b  =  J ic4ia  -  (^  +  p)2  (dropping  term) 

=  i(£  +  p)  J  1  -  ^  -  2 
F  2(£+p)2 


P) 


Then 


tan  ph  M  i 


1  + 


1  -  JBiS. - 

2  (tr  +  P)2 


icqggh 


2  (ir+ph)' 


Let  ph  =  ich,  then 

tan  ph  =  i  tanh  ch 

tanh  ch  =  1  +  ia^g- — ■, 
2  (ir+ich)' 


(A -22) 


(a -23) 
(A -24) 


(A -25) 


(A -26) 


(A -27) 


(A -28) 


-VI  - 


Using  the  relation 

i  „"2ch 

tanh  ch  =  (A-29) 

we  obtain 

e-2ch - J3t2h? -  (A -30) 

4(7r+ich)  +  iajjah* 

To  solve,  it  will  be  necessary  to  apply  an  iteration  procedure, 

2 

using  numerical  values.  At  f  =  1  cps,  icqach*  =  i(0.64). 

Guessing  the  denominator  to  have  the  value  40,  we  have 

-2ch  _  iO.64 

e  = - 2RT" 

2ch  =  log  40  -  log  0.64  +  iTr/2 

c  =  1/h  (2.285  +  i0.785)  =  lA(2.4l  exp  il8.9°)  . 

(A -51a) 

Using  this  value  for  c  in  the  right  hand  side  of  (A-50),  we  ob¬ 
tain  as  the  next  approximation  to  c: 

c  =  1A  (2.1  +  il.56)  =  1/h  (2.61  exp  i56.6°)  (A-51b) 

The  next  iteration  gives 

c  =  lA  (1.88  +  il.71)  =  1A  (2.54  exp  i42 .5°) ,  (a -51c) 

the  next 

c  =  1A  (1.774  +  il.705)  =  1A  (2.46  exp  i45.8°), 

(A-51d) 

and  the  next 

c  =  1A  (1-759  +  il.675)  =  1A  (2.41  exp  i45.9°). 

(A-51e) 


Evidently  we  are  coming  close  to  a  good  approximation,  and  it 
will  suffice  to  use  this  last  figure.  Since  p  =  ic, 
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p  =  i/h  (2.41  exp(i43.9°))  =  1/h  (2 .41  exp  (U33.90) 

=  ttA  (.767  exp  (il33-9°) ) 

and 

Y  =  v  /h  +  p 

Y  =  7T/h[l  +  .767  exp  (1133-9°) 3 

Y  =  7r/h[.468  +  i.552] 

Y  =  IT  Ail  .723  exp  (149.7°)]  - 

Thus  the  magnitude  of  y»  which  at  1  cps  would  be  nearly 
7r/h  for  perfectly  conducting  boundaries,  is  not  greatly  modified 
in  this  case.  The  actual  attenuation  is  reduced  from  ir/h  to 
0.468  vAi,  and  a  phase  shift  o*f  0.552  ir/h  radians  per  meter  is 
introduced.  In  the  A -plane,  this  corresponds  to  the  pole 
moving  from  A^  =  ~  ^  to  ^2  =  E"  (®*552)  -  ^  (0.468). 

7T  2 

At  f  =  1000  cps,  ioyc  »  (j^-  +  p)  .  It  is  easily  found  that 

TT  1 

p  w  - — - 

*  l  +  JI  25-3 
-  (1-i)  9-8  x  10"7 

which  compared  with  2.85  x  10”^  is  fairly  small. 

We  can  summarize  the  above  by  saying  that  if  TE  and  TM 

modes  were  excited  in  the  case  of  perfectly  conducting  boun- 

,  4 

danes  (for  the  frequency  range  1-1000  cps  and  h  =  9  x  10 

meters),  their  amplitudes  would  be  attenuated  with  distance 

much  more  severely  than  the  TEM  mode.  Hence  one  would  expect 

that  they  would  not  be  significant  except  within  short  distances 

from  the  antenna.  It  has  been  shown  that  this  conclusion  is 

still  applicable  when  the  boundaries  are  assigned  a  conductivity 

of  10"^  mhos /meter. 


APPENDIX  B 


The  equations  involved  in  the  calculations  of  the  field 
components  are  presented  in  a  simplified  form  to  facilitate 
their  calculation.  In  these  formulae,  terms  which  are  negli¬ 
gible  compared  with  others,  over  the  frequency  range  1-1000  cps 
have  been  discarded.  The  range  of  validity  of  these  formulae 
is  further  specified  by  parameter  values  listed  under  (2-5-20). 
For  example,  one  cannot  let  h  tend  to  infinity  without  invali¬ 
dating  some  conditions  under  which  the  formulae  are  obtained. 

In  this  way,  it  has  been  found  that 


r(\>) 


ifi  ,k-.h 

.  1  2...  . 

k1(l+ik5h) 


(B-l) 


or  in  terms  of  more  fundamental  parameters 


(B-2) 


Another  quantity  involved  in  the  calculations  is  approximated 
by 


i2X  6.6,h 

(k1  +  Z9)'  as  -  - ,  -  --2  ... 

A.A  CJ1(65-/35h) 


(B-3) 


In  this  formula,  the  second  term  in  the  parenthesis,  in  the 
denominator,  is  negligible  at  lower  frequencies  and  affects 
only  the  second  significant  figure  at  1000  cps.  If  it  is 


dropped,  the  formula  reduces  to 


(k,  +  z2) '  as  - 

X=XQ 


i2X  h 
o 

CDC 


(b-4) 
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Therefore  the  mode  solution  field  components  in  region  1 
arising  from  the  vertical  electric  dipole  can  be  written 
approximately  as 


E 


(v,m) 

lp 


E 


(v,m) 

lz 


H 


(v,m) 

10 


*  Idl 


«  Idl 


w  Idl 


iX  k.coe 
o  i 

4aL2h 


o 


2 

X^  a>e 


o 


4aj^h 


o 


H 


H(2)cXoC)eikl<=-“o) 

H<2>(*0p)eikl<z-z°> 

P,(X0p)eikl(z'zo) 


(B-5) 

(B-6) 


(B-7) 


in  which  M  has  been  replaced  by 
Idl 

4ir  o-^ 


and  Idl  is  taken  to  be  the  moment  of  the  elemental  dipole. 

The'  formulae  for  mode  solutions  in  region  2  are  likewise 
readily  simplified  in  this  way. 

In  the  case  of  the  horizontal  electric  dipole,  the  follow¬ 
ing  simplified  mode  solutions  can  be  written,  using  the  approxi¬ 
mations  previously  indicated: 


E  (h,m) 
Ip 


H 


(h,m) 

10 


-  Idl 


daL(l  +  ik-jh 
r  (^ ) , 


jcos  0[h <2  )  (XQP  )  -Hp  >  (XoP  )]eikl  (z  -zo) 


lXo  k3Hl  (Xop)  ^  ik, (z-z  ) 

m  Idl  - -  cos  0  e  1'  o' 

4a]_k1(l  +  ikjh ) 


(B-8) 

(B-9) 


Idl 


iX  2k,[H^2'(X  p)  -  hI2  ^  (X  p)l  ..  ,  v 

0.  ?L  °  °P/ _ 2  V  °P/j  cos  0  elkl(2-zo) 

8k, (1  +  ik,h) 

*  (B-10) 


-X  - 


In  region  2,  S(X,z),  defined  subsequent  to  equation 
(4.1-11)  can  be  written 

S(X,z)  =  -  xr(x)f1(x,z)  (B-ll) 

where  f-^X/z)  =  j3jCos  02z  +  SjJ^sin  &2Z  * 

Also  r(X,z)  =  -  Xr(X)f2(X,z)  (B-12 ) 

where  f2(X,z)  =  /3-^sin  “  6^2cos  £2Z  * 

The  mode  solution  field  components  in  region  2  can  therefore 
be  approximated  by 


(h,m) 


X  2k,e_lklzo  cos  0 


Idl  *0-3*  M2>(X  p)-h(2)(X  p)L  ,  , 

8a1k1(l  +  ik^h)  L  o  v  op '  2  K  op/Jf1(Xo,z) 

(B-13) 

An3k^e"lklZ°  COS  0  (1\ 


H 


(h,m) 

20 


4o1/B<0>k1(l  +  ik?h) 


ioae  X  2k^e_:''klzo  cos  0 


-  I <31  -  W2)(X  n) 

8a1^2°^k1  (1  +  i*^)  L 

~&2  ^Aop0f2^Xo'z) 


(B-15) 


with  values  of  the  parameters  as  given  in  2.5-20)  it  is  found 
that, as  a  function  of  frequency 


p(°>  *2.34x  IQ"8  PA  e-1  **  (B-l6) 

The  above  approximate  formulae  represent  a  considerable 
simplification  over  the  complete  exact  equations  and  should 
greatly  facilitate  any  computation  which  might  need  to  be  made 
concerning  the  field  components. 
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